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I. INTRODUCTION 



Technological progress in the last two decades has brought about 
the microminiaturization of many electrical components and systems . 

The natural consequence of the transistor was the monolithic integrated 
circuit (IC) . The future promises not only medium-scale circuit inte- 
gration (MSI) but even large-scale integration (LSI) . To make electrical 

filter networks compatible to this trend, the development of inductorless 

< 

filters, better known as active RC filters, came about. 

A. DISADVANTAGES OF INDUCTORS 

Inductance-capacitance (LC) filters were found reliable and easy 
to design at high frequencies but at frequencies below 1000 Hertz, they 
were found limited to a maximum quality factor, Q, of about 100. 

Higher Q's were possible only with large and bulky inductors. Attempts 
to reduce inductor size with high-permeability materials and ferrite 
cores met with failure. Inductor quality deteriorated rapidly at lower 
frequencies . 

There were other disadvantages that came with the use of in- 
ductors. Because of their associated magnetic fields resulting in 
mutual coupling effects and because of their nonlinear behavior, they 
created complications in practical circuit applications. They gave the 
circuit designer additional problems because of incidental dissipation 
resulting from their core loss and winding resistance. Furthermore, 
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inductors were found to be prone to noise pickup. Because of all these 



reasons, the concept of inductorless filters became very attractive for 
many applications. 

B. DEFECTS OF PASSIVE RC FILTERS 

Using resistors and capacitors to meet certain filter specifications 
was found possible but such passive RC filters suffered from two 
defects: 

1. They introduced considerable loss in the passband. 

2. For a given filter specification, they were more complex 
than an equivalent RLC filter. 

C. EVOLUTION OF THE ACTIVE RC FILTER 

It seemed, then, that the answer to those difficulties mentioned 
in the last section was not simply eliminating inductors but replacing 
them with active elements. The basic ideas of associating active 
elements like high- gain and phase-inverting amplifiers with passive 
networks were developed from feedback-amplifier theory [Ref. 1] in 
1945. Guillemin [Ref. 2] in 1949 suggested the use of vacuum tubes 
in conjunction with passive RC elements in a circuit employing the 
principles of feedback. Dietzold [Ref. 3] came up with a patent for 
his RC filter employing a stabilized amplifier in 1951. In 1954 it was 
Linvill [Ref. 4] who finally gave this type of circuit a name - RC 
Active Filters . 
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Since the publication of Linvill’s article, many other synthesis 
methods have appeared and the study of active RC filters has expanded 
considerably. Four active elements have received the most attention: 

1 . Controlled sources 

2. The negative -immitance converter 

3 . The gyrator 

4 . The operational amplifier 

D. ADVANTAGES OF THE OPERATIONAL AMPLIFIER 

There are many reasons for the considerable attention given to 
the operational amplifier. Its biggest advantage, perhaps, over the 
other three active elements is its easy off-the-shelf availability at 
prices that have gone down with the price of the transistor. Because 
of its solid-state construction (transistor and FET differential-input 
configurations) [Ref. 5], it lends itself perfectly to the current trend 
toward microminiaturization. Because of its very high gain (ideally 
infinite), the circuit using it as the active element has very low 
sensitivity to changes in gain. Because of its very low (ideally zero) 
output impedance, cascading low-order functions without the necessity 
of isolating amplifiers is easily accomplished. And, finally, because 
of its very high (ideally infinite) input impedance, resulting in virtual 
ground at its input terminals , circuit analysis and synthesis are both 
simplified . 
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PREVIEW AND PURPOSE 



Chapter III will show the type of analysis involved when oper- 
ational amplifiers are used. The subsequent chapters will show the 
different synthesis techniques that have evolved over the years. These 
techniques have been classified into distinct categories for the purpose 
of unified development and easy comparison. A common starting point 
for all these techniques, namely, the generation of an appropriate 
transfer function to be synthesized, will be developed along the lines 
of modern filter theory in the next chapter. 

It is the purpose of this paper to present a unified and compre- 
hensive study of active RC filters using operational amplifiers. 
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II. MODERN FILTER THEORY 



This chapter will present the development of Modern Filter Theory 
and provide a common starting point for the different synthesis methods - 
the determination of the desired transfer function, H(s) . 

A. HISTORY AND EVOLUTION 

Since the discovery of the electric-wave filter by Campbell and 
Wagner in 1915, filter theory has evolved along two different lines. 

One is known by the name Classical Filter Theory and the other Modern 
Filter Theory. 

The classical theory is more commonly known as the Image- 
Parameter Theory of Filter Design. Originated in the 1920' s by Zobel, 
this theory assumes that, in the design of a filter, its load impedance 
is matched to its image impedance. In practice, however, this 
assumption proves to be inaccurate because most loads are constant- 
value resistances while the image impedance is frequency dependent. 
Consequently, design methods based on this theory involve trial and 
error and often, final adjustments are required to meet design specifi- 
cations . 

Modern filter theory, on the other hand, consists of techniques 

which are more ideal in terms of meeting exact specifications . 

Developed in the 1930's through the efforts of several individuals, 

prominently Norton, Foster, Cauer, Bode, Brune, Guillemin and 

Darlington, this theory has gained considerable acceptance. Essentially, 

19 



it involves the approximation of given specifications with a rational 
transfer function and the realization of this function through the use of 
synthesis techniques. Since such synthesis procedures are analytical 
and exact, no trial and error is involved in designing a filter from its 
transfer function . 

It is along the lines of modern filter theory that the synthesis 
techniques in this paper follow. 

B. TYPES OF FILTERS 

Filters have numerous and varied applications but they all boil 
down to the basic operation of passing through a desired signal and 
rejecting other unwanted signals. 

Every filter is characterized by a range of frequencies over which 
the attenuation of a signal is negligible and a range of frequencies 
over which the attenuation is severe. Such regions of frequency passage 
are termed passbands and the regions of attenuation are called stop- 
bands. In physically realizable filters, there is also a range of 
frequencies between passband and stopband over which the attenuation 
builds up from negligible to severe. Such bands are called transition 
regions. Filters are classified commonly according to their passband 
location in the frequency spectrum . 

The four basic filter types based on this classification are: 

1. Low-Pass Filter: A low-pass filter is a two-port which 

passes all frequencies less than a certain cutoff frequency u) and 

c 

attenuates all frequencies greater than . 
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High-Pass Filter; A high-pass filter is a two-port which 



passes all frequencies greater than a certain cutoff frequency iO and 
attenuates all frequencies less than cj 

c 

3. Bandpass Filter: A bandpass filter is a two-port which 

passes all frequencies between certain lower and upper frequency 

limits, CJ. and co , and attenuates all frequencies outside these 
1 u 

limits . 

4. Band-Reject Filter: A band-reject filter is a two-port which 

attenuates all frequencies between certain lower and upper frequency 
limits, and passes all frequencies outside these limits. 

The effect of a filter on the passage of frequencies is presented 
graphically in terms of gain or magnitude and phase characteristics. 

In their ideal form, the magnitude and phase characteristics of the 
four basic types are shown in Figure 2-1. H(s) is the filter transfer 
function representing the ratio of filter output to input transforms, 
where s is the complex frequency variable. In their more practical 
form, the attenuation characteristic of the four basic types are shown 
in Figure 2-2 . 

C. THE FILTER DESIGN PROBLEM 

The ideal filter is one which passes all frequencies without 
attenuation within its passband, completely attenuates all frequencies 
within its stopband and has a linear phase response in its passband. 

Such ideal filter characteristics , however, are physically un- 
realizable. No quotient of rational, finite-degree polynomials can 
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Figure 2-1. Ideal magnitude and phase characteristics of 
(a) a low-pass filter (b) a high-pass filter 

(c) a bandpass filter (d) a band-reject filter 





Figure 2-2. Attenuation characteristics of (a) a low-pass 
filter (b) a high-pass filter (c) a bandpass filter 
(d) a band-reject filter 
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represent exactly the ideal magnitude and phase characteristics. In 
fact, the desirable abrupt transition from zero to infinite attenuation 
with negligible change in frequency in the gain characteristic would 
wreck havoc to a filter's transient response because of the infinite 
phase shifts and infinite lags in frequencies within the passband, or 
in effect, a drastic departure from the desired linear phase characteristic. 
On the other hand, a more gradual fall off of magnitude with frequency 
change results in a more acceptable phase response. In general, then, 
the fundamental problem in modern filter design is in finding the transfer 
function H(s) that best approaches the ideal magnitude and phase 
characteristics. For other particular filter applications, it can happen 
that one characteristic's requirements are relaxed in favor of more 
stringent requirements in the other. This process of finding the best 
H(s) depending on the desired filter application is known as the approx- 
imation problem. Once this H(s) is known, the other half of the modern 
filter theory problem is the exact realization of the filter network using 
synthesis techniques . 

D. FILTER CHARACTERISTICS 

In this section, three approximations to the ideal low-pass filter 
will be discussed briefly to give the reader an idea of the kind of 
compromise involved in meeting both frequency- and time-response 
specifications. These are the Butterworth, Chebyshev, and Bessel or 
Thompson filters . Each approximates the characteristics of the ideal 

filter in a different way and each has its area of application. 
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There are other types of filters such as the Cauer or elliptic- 



function filters, the monotonic Legendre filters and the filters with 
parabolic pole distribution, more commonly known as the P filters. 
Detailed information about their frequency and time responses as 
compared with those of the types to be discussed in the following 
paragraphs may be found in Refs. 6, 7, and 8. 

Only the normalized low-pass version of each type will be shown 
since the high-pass, bandpass and band-reject transfer functions 
exhibiting the same characteristics of the low-pass prototype can be 
easily derived using the appropriate frequency transformations . Such 
transformations will be discussed in the next section. 

1 . Butter worth Class 

Butterworth filters are characterized by the normalized 
magnitude-squared function: 



H(j ^ ) 



1 + 






2n 



( 2 - 1 ) 



where n is the number of poles in the network, is the frequency of 

interest and CO is the cutoff frequency. ^ = 1 for the normalized 

c c 

case. From equation (2-1), it can be observed that regardless of the 
number of poles, the output at the cutoff frequency is always 3 dB 
below that at CiJ = 0. The derivative of this function shows that the 
slope of the response curve is zero at = 0. For this reason the 
Butterworth filter is called a maximally flat amplitude filter. At 



cutoff {CJ = CtJ ) , the slope of the response curve is (-n/2), showing 
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an increasing rate of attenuation with increasing number of poles. 
Beyond cutoff the filter response falls at an approximately constant 
2 On dB / decade . 

By using the stable poles of the magnitude-squared 
function, it can be shown that the Butterworth low-pass transfer 
function is of the form: 

H(s) = K/Bn(s) (2-2) 

where K is a constant multiplier and the polynomials Bn(s) are called 
the Butterworth polynomials. They are given in Table 2-1 up to order 5. 
Their roots, all of which lie on the unit circle, are listed in Table 2-2. 
The low-pass amplitude response for different values of n is shown in 
Figure 2-3. It can be noted that the amplitude response flatness over 
the entire passband increases with the number of poles. 

The phase characteristic, however, is not very linear. 
Consequently, the time delay changes with change in frequency. The 
effect of this is an overshoot and some ringing in the filter's response 
to a step function. The amount of overshoot increases with the number 
of poles. The rise time of the output is satisfactory but increases 
rather slowly with increasing n. 

2 . Chebyshev Class 

The Chebyshev filter is characterized by the magnitude- 
squared function: 





H(j^) 



(2-3) 







/ 



m 



T 



A, 



n 

1 

2 

3 

4 

5 



Table 2-1. The Butterworth Polynomials 



s +1 

+ '/2 s + 1 

+ 2s^ + 2s + 1 = (s + l)(s^ + s + 1) 

s^ + 2.613s^+3.414s^+2.613s+l = (s^+0 . 765s+l) (s^+1 . 848s+l) 

5 4 3 2 

s + 3.236s +5. 236s +5. 236s +3.236s+l 

2 2 

= (s+1) (s +0.618s+l)(s +1.618s+l) 



Table 2-2. Roots of the Butterworth Polynomials 



n = 1 n = 2 


n = 3 


n = 4 


n = 5 


-1.0000000 -0.7071068 

+ j0. 7071068 


-1.0000000 


-0.3826834 
+ jO. 9238795 


-1.0000000 




-0.5000000 
+ j0. 8660254 


-0.9238795 
+ jO. 3826834 


-0.3090170 
+ j0. 9510565 








-0.8090170 
+ jO. 5877852 
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FREQUENCY in radian/second 
Figure 2-3. Butterworth low-pass amplitude responses for 
various values of n. 
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where 6 is the ripple parameter and Tn(u>/ uj ) is the Chebyshev 

c 

Polynomial of order (and degree) n. 

The transfer function can be shown to be: 

H(s) = K/Vn(s) (2-4) 

where K again is a constant and Vn(s) is formed from the left half 
plane zeros of the denominator of equation (2-3) . Each value of £ 
yields a different set of Vn polynomials . The coefficients of these 
polynomials corresponding to two values of £ up to n = 5 are given in 
Table 2-3; their roots are given in Table 2-4. 

The magnitude response of this class of filters is dis- 
tinguished by an equal-magnitude ripple in the passband and a maximum 
rate of fall off beyond cutoff. The ripple height, ^ , is given by: 



6 = 1 - 






(2-5) 



The sum of the number of maxima and minima in the passband is equal 
to the order of the filter. A fifth-order Chebyshev low— pass filter 
response is shown in Figure 2-4. The response at the cutoff frequency 
is always that of a minimum. Thus, a Chebyshev filter with + 1-dB 
ripples will be - 1 dB down at cutoff. For the + 3-dB ripple case, both 
Butterworth and Chebyshev filters will have the same cutoff frequency. 
In general, the 3-dB-down point will be at: 

cosh cosh "1/e Cl) 



CJ 



3dB 



n 



(2-6) 
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n 

1 

2 

3 

4 

5 

n 

1 

2 

3 

4 

5 



Table 2-3. Polynomials for Chebyshev Filters 



(a) 1/2 dB Ripple 



s +2.863 

s^ + 1.425s +1.516 

s^ + 1.253s^ + 1.535s + 0.716 = (s+0 . 62 6) (s^+0 . 626s+l . 142) 
s"^ + 1.197s^+1.717s^+1.025s+0.379 = (s^+0 .351s+l . 064) 

(s +0 . 84 5s+0 . 3 56) 

s^ + 1.172s^+1.937s^+1.309s^+0.753s+0.179 

= (s+0 .362) (s^+0 . 224S+1 . 036) (s^+0 . 586s+0 .477) 

(b) 1 dB Ripple 

s + 1 . 965 

s^ + 1.098s + 1.103 

+ 0.988s^ + 1.238s + 0.491 = (s+0 . 494) (s^+0 . 494s+0 . 994) 
+ 0.953s^+1.454s^+0.743s+0.276 = (s^+0 . 279s+0 . 987) 

(s^+0.674s+0.279) 

s^ + 0.937s'^+1.689s^+0.974s^+0.581s+0.123 

= (s+0.289)(s^+0. 179s+0.988)(s^+0.468s+0.429) 
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Table 2-4. Roots of Chebyshev -derived Polynomials 
(a) 1/2 dB Ripple (C = 0.3493114) 



-2.8627752 -0.7128122 


-0.6264565 


-0.1753531 


-0.3623196 


+ jl. 0040425 




+ jl. 0162529 






-0.3132282 


-0.4233398 


-0.1119629 




+ jl. 0219275 


+ jO. 4209457 


+ jl .0115574 








-0.2931227 
+ j0. 6251768 



(b) 1 dB Ripple (€ = 0.5088471) 



-1.9652267 -0.5488672 


-0.4941706 


-0.1395360 


-0.2894933 


+ j0. 8951286 




+ jO. 9833792 






-0.2470853 


-0.3368697 


-0.0894584 




+ jO. 9659987 


+ jO. 4073290 


+ j0. 9901071 








-0.2342050 








+ j0. 6119198 
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MAGNITUDE 




FREQUENCY in radian/second 
Figure 2 - 4 . A fifth-order Chebyshev approximation to a 
low-pass filter. 
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The rate of attenuation beyond cutoff is initially greater than 20n dB/ 



decade and increases with the magnitude of the allowed ripple and the 
number of poles . 

The phase response is even more nonlinear than that of the 
Butterworth class and this nonlinearity increases with the allowed 
ripple and the number of poles . 

3 . Bessel Class 

The Bessel filter exhibits maximally flat time delay 
(linear phase) and subsequently its response to a step function shows 

I 

less than one percent overshoot. 

The amplitude response, on the other hand, is not flat 

but decreases monotonically in the passband. The cutoff frequency, 

U) / is not the 3-dB-down point but is defined in terms of the zero- 
c 

frequency time delay, t^, as: 



CJ = 
c 



(2-7) 



The location of the half-power point, OJ is a function of the 

order and for n > 3 is given approximately by: 



s/o.69315(2n-l) 



<^3dB ^ VU.by3i5izn-i; (2-8) 

The Bessel class of filters is represented by the low-pass 



transfer function: 

H(s) = K/Hn(s) 



(2-9) 
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where K is a constant and Hn(s) are the Bessel Polynomials , Their 
coefficients up to n = 5 are given in Table 2-5; their roots are given 
in Table 2-6 . 

4 . Comparison of Classes 

Figure 2-5 contrasts the magnitude characteristics of the 
three filter classes. For a given pole order, magnitude flatness 
changes drastically as one goes from Butterworth to Bessel to 
Chebyshev. Sharpness of cutoff improves from Bessel to Butterworth 
to Chebyshev. Linearity of phase deteriorates rapidly in the same 
order. Consequently, percentage of overshoot of step input response 
increases while rise time improves in this same order. A comparison 
of these last two figures of merit is given in Table 2-7. From this 
comparison of approximations, it can be seen dramatically how the 
choice of filter type must depend on the relative importance of ampli- 
tude and phase responses to the particular application desired. 

E. FREQUENCY TRANSFORMATIONS 

High-pass, bandpass and band-reject filters exhibiting the same 
characteristics as their low-pass prototype may all be obtained by 
making the following frequency transformations; 

Low-pass to high-pass: 

(jL^ C 

Sn = ^ (2-lOa) 

s 
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Table 2-5. The Bessel Polynomials 

p (s) = 1 
o 

p^(s) = s + 1 
2 

p^ (s) = s + 3s + 3 

p (s) = s^+6s^+15s+15 = (s+2 .322)(s^+3.678s+6.460) 

O 

p^(s) = s^+10s^+45s^ + 105s+105 = (s^+5 . 792s+9 . 140) (s^+4 . 208s+l 1 . 488) 

p (s) = s^+15s +105s^+420s +945s+945 
o 

= (s+3.647)(s^+6.704s+14.272)(s^+4.649s+18.156) 

p (s) = (2n - 1) P„ 1 (s) + s^p _(s) 
n n-i n-z 

Table 2-6. Roots of the Bessel Polynomials 
n 

1 - 1.0000000 

2 - 1.5000000 + jO. 8660254 

3 - 2.3221854; - 1.8389073 + jl. 7543810 

4 - 2.8962106 + jO. 8672341; - 2 . 1037894 j^j2. 6574180 

5 - 3.6467386; - 3.3519564 + jl . 7426614; - 2 .3246743 + j3 .5710229 
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fifth-order filters . 
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N 



I 

I 



Table 2-7 



% Overshoot 


n 


Butterworth 


Bessel 


Chebyshev 1/2 dB 
Ripple 


Chebyshev 1 dB 
Ripple 


1 


0. 


0. 


0. 


0. 


2 


4.3 


0.43 


10.5 


14.6 


3 


8.15 


0.75 


9.2 


6.6 


4 


10.9 


0.83 


17.8 


22.2 


5 


12.8 


0.77 


13.1 


10.2 



Rise Time (10 to 90%) 



n 


Butterworth 


Bessel 


Chebyshev 1/2 dB 
Ripple 


Chebyshev 1 dB 
Ripple 


1 


2.2 


2.2 


2.2 


2.2 


2 


2.15 


2.73 


1.85 


2.0 


3 


2.29 


3.07 


1.81 


1.97 


4 


2.43 


3.36 


2.0 


2.03 


5 


2.56 


3.58 


1.98 


1.97 
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Low-pass to bandpass: 



Sn = 



C(J o 

BW 



Oj o 



COc 



(2-lOb) 



Low-pass to band-reject: 



Sn = 



BW 



t0o(- 



(J 









(2-lOc) 



where S represents the normalized low-pass complex frequency 
n 

variable, QJ ^ the center frequency and BW the bandwidth of the 



bandpass and band-reject filters. Center frequency is defined as: 




while bandwidth is 

BW = (2-12) 

where is the lower frequency limit and the upper frequency 
limit of a bandpass or band-reject filter. 

Strictly, bandwidth is defined for only the low-pass and band- 
pass filters. For the low-pass one, the bandwidth is equal to cO • 

c 

For the bandpass case, it is represented by equation (2-12). Normally, 
the half-power points in the amplitude response are used to indicate 
these lower and upper frequency limits. In this case, the bandwidth 
is called the half-power bandwidth. 

These transformations are performed on the normalized low-pass 
transfer function H(s) . It can be noted that the center frequency is 
the geometric mean of the upper and lower frequency limits and will 
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not be necessarily in the center of the bandwidth, BW. Consequently, 



final designs may be expected to be slightly asymmetric with respect 
to the intended center frequency. 

F . FREQUENCY SCALING 

It can be noted that the low-pass transfer function, H(s), of the 

different filter classes has been presented on a normalized frequency 

basis (i.e., r.) =1 rad / sec). One reason for this is the ease of 

^c 

working with simple normalized coefficients. If it is desired to shift 
the response of the filter up to a desired cutoff frequency, CU^' 
frequency scaling is easily accomplished using the relation: 

Sn = — 7 — (2-13) 

LJc 

Another reason is flexibility. In frequency transformations, the equations 
have already been conveniently scaled to operate on the normalized 
low-pass transfer function. For these reasons, filter design hand- 
books and tables [Refs. 9, 10] publish only the poles and zeros of the 
normalized low-pass transfer function of the different classes of 
filters for different orders . A frequency transformation procedure 
operating directly on poles and zeros given in Ref. 9 was shown by 
J. Tow [Ref. 11] and is contained in Appendix A. 

G. THE FILTER DESIGN PROCEDURE 

At this point, a filter design procedure may now be outlined. It 
will consist of two phases: 
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1 . 



Determining the transfer function, H(s), from the given 



specifications . 

2 . Realizing the transfer function by means of available 
synthesis procedures . 

Given a set of specifications consisting of a desired amplitude 

response and a desired phase response or a desired time response to a 

step input, the best approximation from the different filter classes is 

selected. If the filter desired is low-pass, then all that is involved 

is frequency scaling to the desired cutoff frequency, ^ . If the 

c 

filter desired is high-pass, bandpass, or band-reject, then the 
appropriate frequency transformation is done to come up with the re- 
quired transfer function, H(s). 

Once this transfer function has been decided upon and determined, 
the active filter network ma’y be synthesized using any of the different 
synthesis procedures that will be discussed in the following chapters. 

It is useful, at this point, to observe that if the general second- 
order transfer function 



H(s)=^ 




(2-H) 



s + b, s + b 
1 o 



can be written as 



H(s) = K . 




o 



n 



2 



2 



(2-15) 
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where K is an arbitrary gain factor, is the resonant frequency, Q 

the selectivity factor, the "notch frequency," and Q the "notch 

n n 

attenuation factor," then the four different filter types may be con- 
veniently written as follows: 

Low-pass; 



H(s) = 



+ ( coyQ) s + 



High-pass: 



(2-16) 



H(s) = 



K.S 

+ { Cj /Q) s + U) ^ 
o o 



(2-17) 



Bandpass: 



H(s) 



s 

o 

s^ + ( U)/Q) s + 

’ o o 



(2-18) 



Band-reject; 



H(s) = K . -2 ^2 

s + ( CL>yQ)5 + 

where, in each case, the numerator polynomial of equation (2-15) has 
been modified to conform with the transfer function characteristics of 
each filter type . 

The second-order transfer function is emphasized because not 
only does it serve as a simple model to realize using the different 
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synthesis techniques but, more important, the current trend in realizing 
filter transfer functions of any desired order is to break them up into 
first- and second-order product terms , realize each of these terms 
individually using first- and second-order building-block networks 
and cascade these to produce the desired transfer function. This con- 
cept will be treated in greater detail in Chapter VII of this paper. 



H . EXAMPLES 

Several examples to illustrate the first phase of the design pro- 
cedure outlined in the preceding section follow. 

Example 2-1. A second-order Butterworth filter operating at a 
low-pass cutoff frequency of 10 Hz is to be realized. Specify its 
transfer function. 

The desired transfer function can be specified directly by using 
equation (2-16): 



H(s) = 



3947.84 



s + 88.84s + 3947.84 



( 2 - 20 ) 



whereK=l,6; = 2 TT (10) and l/Q = 1 .414 . 
o 

Computer Program 1 contained in Appendix B was used to give a 
magnitude-frequency plot of this low-pass filter. The plot is shown 
in Figure 2-6 . 

Example 2-2. Specify a transfer function for a low-pass filter 

with a cutoff frequency of 10 Hz. Its magnitude response is to be 

essentially flat in the passband. A fall-off rate of 100 dB / decade 

beyond cutoff is desired. Linearity of phase is desirable but not necessary. 
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FREQUENCY in Hertz 

Figure 2-6. Magnitude vs Frequency Plot of Filter Transfer 
Function of Example 2-1. 
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A fifth-order Butterworth filter is chosen. Its normalized low- 
pass transfer function from Table 2-1 is 



H(s) = -z i 7 (2-21) 

s + 3.236s + 5.236s + 5.236s + 3.236s + 1 

or, in factored form. 



H(s) = — — • -r • T (2-22) 

s+1 s + 0.618s +1 s + 1.618s + 1 

Frequency scaling is called for and the substitution 

Sn = s/2 IT (10) (2-23) 

is applied to equations (2-21) and (2-22) resulting in the desired 
low-pass filter transfer function 



H(s) = 



979261972.65 

s^ + 203.328"^ + 20670. 89s^ + 1298790. lls^ 

+ 50434489.00s + 979261972.65 



(2-24) 



or, in factored form. 



H(s) = 



3947.84 



3947.84 



62.83 . 

s+62.83 s^+38.83s+3947.84 s^+10 1 . 66s+3947 . 84 



(2-25) 



Its magnitude response plot is shown in Figure 2-7. 

Example 2-3. A high-pass filter with a 3-dB ripple in the pass- 
band operating at a cutoff frequency of 10 Hz is desired. Phase response 
is not critical. Specify its transfer function. 

A 3-dB ripple Chebyshev filter would satisfy the requirements. 

Its normalized low-pass transfer function [Ref. 10] is 
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Figure 2-7. Magnitude vs Frequency Plot of Filter Transfer 
Function of-Exarriple' 2 -2 . 
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H(s) = 



0.7079478 



(2-26) 



+ 0.6448996s + 0.7079478 
The frequency transformation relation 

Sn = 2 TT (10)/s (2-27) 

is applied to equation (2-26). The resulting frequency-transformed 
transfer function is 



H(s) = 



(2-28) 



s + 57.24s + 5576.46 



The magnitude-frequency plot using Computer Program 1 is shown in 
Figure 2-8. 

Example 2-4. Specify a second-order bandpass transfer function 
with center frequency at 10 Hz and a Q of 50. 



Application of equation (2-18) gives 

LX)0 - 

“S' 

62^3s , 



i.. 






1- 









H(s) = 



s^ + (62.83/50)s + (62.83)^ 



or 

H(s) = — r 62^-83 s . ^2-29) 

s + 1.257s + 3947.84 

The magnitude response plot is shown in Figure 2-9. 

Example 2-5. A bandpass filter is desired with center frequency 
at 10 Hz and Q of 100. Essentially linear phase is desired within the 
passband. Specify its transfer function. 
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Figure 2-8 • Magnitude vs Frequency Plot of Filter Transfer 
Function. ofiiE xam ple:2- 3 . 
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Figure 2-9. Magnitude vs Frequency Plot of Filter Transfer 
Function of Example 2-4. 
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The linear phase requirement calls for a Bessel response, so the 



appropriate frequency transformation 



Sn = 



1 



100 



2 TT (10) 



2 Tf (10) 

s 



(2-30) 



is applied to the normalized low-pass Bessel transfer function from 
Table 2-5, 



H(s) = 



s + 3s + 3 



(2-31) 



A fourth-order bandpass transfer function results: 



H(s) = ~2 ^ ' (2-32) 

s +1.88S +7896. 86s +744 . 50s+15585441 . 60 

The fourth-order denominator of the transfer function in equation 

(2-32) may be factored ^ into two second-order polynomials to give 



H(s) = ^ • (2-33) 

(s +0.94S+3913.80) (s +0 . 95s+3982 . 18) 

The magnitude-frequency plot of equation (2-32) is shown in Figure 

2 - 10 . 

Example 2-6. A filter to reject signals at frequencies between 9 
to 11 Hz is desired. Second-order Chebyshev characteristics with l/2-dB 
ripples are specified. Generate the transfer function to be synthesized. 



^Subroutine POLRT was used to determine the roots of the fourth- 
order denominator polynomial. 
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FREQUENCY in Hertz 

Figure 2-10. Magnitude vs Frequency Plot of Filter Transfer 
Function of Example 2-5. 
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Using equation (2-11), the appr opriate center fr equency, is 

calculated: oao ^ 

6o)q= / 2 TT (9 .11) = /i98TT . (2-34) 

Bandwidth, BW, is determined by equation (2-12): 

BW = 2 TT (11-9) = 4 TT . (2-35) 

Low-pass-to-band-reject frequency transformation given by equation 
(2-lOc) is then applied to the normalized second-order low-pass 
Chebyshev transfer function from Table 2-3a: 

H(s) = • (2-36) 

s + 1.425s + 1.516 

The resulting transfer function is 



H(s) = 



s"^ + 7816s^ + 152 75291. 3 1 



s'^ + 11.81s^ + 7920. 89s^ + 46165.78s 



(2-3 7) 



+ 15275291.31 



Factored for easy second -order realization, the transfer function 
becomes: 



s^+3908.36 . s^+3908.36 

2 2 
s +5.51S+3420.98 s +6 . 30s+4465 . 19 



(2-38) 



The magnitude response plot of equation (2-37) is shown in Figure 2-11. 
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Figure 2-11. Magnitude vs Frequency Plot of Filter Transfer 
Function of Example 2-6. 
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III. THE OPERATIONAL AMPLIFIER 



Having established a common starting point for the different 
synthesis techniques, it is now important to take a closer look at the 
common active element that is used in all of these methods and the type 
of analysis involved. 

This chapter defines ^ the operational amplifier and shows how its 
distinctive characteristics simplify analysis of circuits using it. These 
same characteristics will be used to justify some fundamental assumptions 
in the synthesis techniques . 

A. DEFINITION 

The operational amplifier is a high-gain, DC -coupled amplifier 
having either a differential or single-ended input. The output is usually 
single-ended with respect to ground. For the differential-input type, a 
circuit model is shown in Figure 3-1. Signals applied to the terminal 
marked + In are amplified by a positive non-inverting gain, + A. Signals 
applied to the terminal marked - In are amplified by a negative - A. The 
output is given by E^= A(E 2 - E^^) . The single-ended amplifier may be 
treated as the special case where + In is grounded. 

The idealized characteristics of the operational amplifier are: 

1 . Gain = CO (a CO ) 

^Courtesy of the Burr-Brown Research Corporation from the "Burr- 
Brown Handbook and Catalog of Operational Amplifiers," 1969. 
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Figure 3-1. Circuit model for the differential-input type of 
operational amplifier 




Figure 3-2. Idealized model of the differential-input type of 
operational amplifier 
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2 . 



E = O when E = E„ 
o 12 

3. Input Impedance = CO (Z^ ^CD ) 

4. Output Impedance = O (Z^ O) 

5. Bandwidth = CO (response time = O) 

When these idealized characteristics are incorporated, the circuit 
model of the operational amplifier reduces to that of Figure 3-2 . This 
idealized model will be used in developing equations for many of the 
basic feedback circuits in later chapters. 

The most common circuit symbols for the operational amplifier are 
shown in Figure 3-3 . 

B. ANALYSIS 

To illustrate the kind of circuit analysis that follows from the 
idealized model, a fundamental inverting circuit such as Figure 3-4 
will be used as an example . ' 

The common feature of these circuits is that the non-inverting 
input is connected to ground. Feedback and input networks are attached 
to the inverting input terminal. In analyzing the circuit, a gain A is 
assumed. This gain is subsequently allowed to approach infinity. 

Since amplifier input impedance is infinite, zero current flows into the 
amplifier; consequently and carry equal currents. This may be 
stated as Ij^ = I^ or 
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(a) (b) 



Figure 3-3. Common circuit symbols for the operational amplifier 
(a) Differential-input (b) Single-ended 




Figure 3-4. Fundamental inverting amplifier circuit 
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The amplifier gain enforces the condition 

E = -AE (3-2) 

o s 

Solving equation (3-2) for E and substituting into equation (3-1) yields 

s 



E 

o 



E^ + A 





(3-3) 



If the assumption A — *• OO is now applied, equation (3-3) becomes 




The overall (closed loop) gain, or transfer function, is 



(3-4) 




Note that the gain includes a sign inversion . This sign inversion will 

be appropriately indicated in the different transfer functions that were 

generated in Chapter II when these are used as examples for realization 

in the following chapters. Note also that its magnitude is determined 

solely by the ratio of external resistors. Another important point is that 

the summing point voltage, E , approaches zero as the gain. A, approaches 

s 

inf in ity . 



E 

E^=- ^ 0;A— ’-CP • (3-6) 

s A 

This condition is described by referring to the summing point as a 
"virtual ground" (since it is held at virtually zero or ground potential) . 
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With the summing point at ground potential, the current through is 

and therefore entirely independent of R^ • Since no current 

flows into the amplifier, the input circuit may be thought of as a source 

of current which must flow through the feedback impedance, in this case 

R . Since one end of R is at ground potential (the summing point), the 

other end must be at a voltage of - I, R_ = E . 

i I o 

The feedback element need not be a resistor but can be any imped- 
ance . This will be seen in the next chapter. In fact, it need not be 
a linear element but can be any element, or set of elements, for which 
there is a linear relation between short-circuit current and terminal 
voltage. Similarly the input circuit is subject to equal freedom. 

Figure 3-5 illustrates a complex network of this type. 

If neither input to the operational amplifier is grounded, the same 
analysis still holds. Consider the network in Figure 3-6. The output 

voltage, E , is 
o 

As A approaches infinity. 

Observe that the operational amplifier's infinite gain forces its two 
inputs to be equal. It will also be noted that since no current flows 
into the operational amplifier 
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Figure 3-5. Complex linear feedback circuit 



R, 



R, 




Figure 3-6. Feedback circuit with neither input to the operational 
amplifier grounded 
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and 



I 



1 






(3-10) 



Simplifying equation (3-10), the voltage transfer function of the network 
is 




^ '"'^2 

R, 



(3-11) 



Regardless of the complexity allowed in the input and feedback 
networks, the same summing point restraints will always be found to 
hold: 



1. The summing point is a virtual ground. 

2 . No current flows into the amplifier; current flowing into the 
input circuit must flow into the feedback circuit. 



C. NATHAN'S MATRIX ANALYSIS 

Since matrix analysis lends itself to more complicated circuits, 
it would also be convenient to know a matrix analysis of circuits con- 
strained by an operational amplifier proposed by Nathan [Ref. 12]. Such 
analysis will be illustrated in later chapters and a brief discussion of it 
is contained in Appendix C of this paper. 
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IV. INFINITE -GAIN SINGLE-FEEDBACK REALIZATIONS 



This chapter presents the first of five categories into which 
synthesis techniques for active RC filters using operational amplifiers 
have been classified. The synthesis procedures involve one infinite-gain 
operational amplifier with a single feedback loop and consequently, the 
significant contributions falling under this general category are presented 
and are collectively known as Infinite-Gain Single-Feedback Realizations. 



A. THEORETICAL DEVELOPMENT 



The basic active filter configuration using one operational ampli- 
fier is shown in Figure 4-1. It employs one RC feedback network, N , 

B 

across the operational amplifier, an RC network, N , in series with the 
inverting input source and another RC network, N^, in series with the 
non-inverting input source . 

The following two-port equations may be written for this active 
configuration: 



^llA^l 


+ 


^12A^3 


= ^21A^1 


+ 


^22A^3 


^ ^11C^2 


+ 


^12C^4 


“ ^21C^2 


+ 


^22C^4 


" ^IIB^S 


-f- 


^12B^6 


“ ^21B^5 




y E 
^22B 6 
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(4-la) 

(4-lb) 



(4-2a) 
(4 -2 b) 



(4-3a) 



(4-3b) 




Figure 4-1. Basic single-feedback active filter configuration 




Figure 4-2 . Differential-input operational amplifier using resistors 
for input and feedback elements 
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Application of one of the summing point restraints developed in 
Chapter III, namely, that no current flows into the operational amplifier 
because of its infinite input impedance, results in: 



1=0, (4-4) 

o 

= 0 . (4-5) 



Because of equation (4-4), I = - I , so that from equations (4-lb) and 

o 0 

(4-3a): 



^21A^1 



V E = - 
^22A 3 






" ^12B^6‘ 



(4-6) 



As a consequence of equation (4-5), equation (4-2b) becomes: 



E 



4 



l21C 

^22C 




(4-7) 



Use of the other summing point restraint, namely, that the summing 
point is a "virtual ground" yields: 

Applying equation (4-8) to equation (4-7) gives: 

^5 " ~ ^21C ^3 ' (4-9) 

^22 C 

^3=ll2IC S 

^22C 

Substituting equations (4-9) and (4-10) into equation (4-6) gives an 

expression for E_: 

6 
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y + V 
IIB ^2 2 A 



(4-11) 



E =E = — V E +/y 
6 o ^21A 1 / ^2 1C 






12B 



22C 



12B 



Note that E^is expressed as a function of the two input voltages and 
the admittance parameters of the three passive networks, N^, and 
N^. To verify this result, note that if the three passive networks were 
each single resistors as shown in Figure 4-2, such that 



^ 21 A " ^^^1 



'^22k ^'^^l 



^ 12 B = - ^/^2 



^ llB =^/^2 



(4-12) 



''21c ~ ^^^3 



^22C ^^^3 ' 



application of equation (4-11) gives: 



R, \ R, / 



(4-13) 



Note that if the non-inverting input, E = 0, equation (4-13) corresponds 
to equation (3-5) and if E^ =0, equation (4-13) corresponds to equation 
(3-11). 

The general expression given by equation (4-11) can now be 
applied to a more restricted case where the non-inverting input voltage, 
£2 = 0 , giving: 



21A 

^12B 



(4-14) 



which is the open-circuit voltage transfer function of the more common 
network configuration shown in Figure 4-3. This is the basic network 
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Figure 4-3. Single-ended operational amplifier in basic single- 
feedback configuration 



Output Network 




Figure 4-4. Active Network configuration of Bradley and McCoy 
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configuration to produce second-order transfer functions. Its importance 



lies in the fact that higher orders are achieved by simply cascading such 
configurations . 

Equation (4-14) is the starting point for all the synthesis methods 
that fall under this category. A more rigorous derivation of this equation 
is contained in Ref. 22. 

B. BRADLEY AND McCOY [1952] 

Writing about driftless DC amplifiers for a Reeves model A- 105 

computer, Bradley and McCoy [Ref. 13] proposed the use of DC amplifiers 

in combination with passive networks to generate complex transfer 

functions as a generalization of summation or integration. Calling Z. 

and Z the short-circuit transfer impedances of the input and feedback 
o 

networks respectively, as shown in Figure 4-4, they used as a basis 
for their work the voltage transfer function, 

^o , (4-15) 




where Z = l/y,o,^ Z. = l/y„,,v in relation to equation (4-14). 
o i z B 1 z lA 

They generated an extensive table of transfer impedance functions 
with their corresponding RC networks and accompanying relations to 
determine element values for the networks. 

Synthesis of a desired transfer function involves choosing appro- 
priate input and output networks such that their transfer impedances , 
when applied to equation (4-15), produce the desired transfer function. 
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Element values of each network may then be chosen by coefficient- 
matching . 

An example of this procedure follows. 

Example 4-1. Realize the low-pass Butterworth filter transfer 
function of Example 2-1, where the negative sign is introduced because 
of operational amplifier inverting action: 



H(s) = 



V^(s) 



3947.84 



s + 88.84 s + 3947.84 



(4-16) 



or, in a more convenient form for this procedure. 



H(s) = - 



(4-17) 



1 + 88 . 84 s + 
3947.84 



3947.84 



Two transfer impedance functions are now taken off the transfer 
impedance table in Ref. 13: * 



(1 + sT) , 



(4-18a) 



where 



A^=R^ 



T = R^C , 



(4 -18b) 



and 



A 



2 



1 + sT, 



1 + sT^ + s 




(4-19a) 
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where 






Ij = 2R^C^ 



= R3C/2 



(4-19b) 



The RC configuration corresponding to equation (4-18) is shown in 

Figure 4-5a; that for equation (4-19) is shown in Figure 4-5b. 

It can be observed that if Z, = A, (1 + sT) and Z = A_ (1 + sT„)/ 

1 1 o 2 z 

(1 + sT^ + s^T^T^), 



H(s) = 



A 2 (1 + sT2)/(1 + sT^ + s T^T^) 
A^ (1 + sT) 



(4-20) 



Setting T = gives 



H(s)= . 

2 

1 + sT^ + s 

Equating coefficients of equations (4-21) and (4-16) yields 



^ 88.84 

1 3947.84 

T T = 

12 3947.84 

Simplified, the relations in equation (4-22) become 



(4-21) 



(4-22) 
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C = T/A 



1 
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2 ' 2 
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Figure 4-5. Transfer impedance networks from the table of Bradley 
and McCoy 
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Figure 4-7. Active network configuration of Mathews and Seifert 
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(4-23) 



^2 ^1 



= 2.245 X 10 



-2 



= 1.127 X lO"^ = T . 

Applying the relations in Figures 4-5a and 4 -5b and arbitrarily choosing 
= 10 K give the component values for the two networks: 

^2 " = 5 KQ 



C= ^- . 127 X 10 . ^ 1.127 juF 

10 ^ 






1 



10 



C = ^- . 245 X 10 . ^ 2.245 jUF 

^ 10 



(4-24) 



and the resulting circuit realization is shown in Figure 4-6, 



C. MATHEWS AND SEIFERT [1955] 

Mathews and Seifert presented a paper [Ref. 14] which described 
three systematic synthesis methods for realizing complex transfer 
functions. The first of these belongs to the single-feedback realization 
category of this chapter. 

Mathews and Seifert approached the synthesis problem in two ways. 
Using the basic feedback configuration shown in Figure 4-7, they 

tt 

assumed an infinite-gain amplifier and initially treated Y (s) and Y„(s) 

A D 

as two-terminal RC driving-point admittance networks to come up with 
the open-circuit voltage transfer function: 
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(4-2 5) 



E (s) Y (s) 

O _ A 

E^(s) ■ Yg(s) 

Consequently, the poles and zeros of each Y (s) and Y (s) must alternate 

A D 

along the negative real axis of the complex-frequency plane and the lowest 
critical frequency of each must be a zero. In effect, the poles and zeros 
of the transfer function must all lie on the negative real axis. However, 
since the transfer function poles are the poles of Y (s) and the zeros of 
Y (s) and since the transfer function zeros are the zeros of Y (s) and the 

D ^ 

poles of Y (s) , two transfer function poles or two zeros may occur to- 
B 

gether and the transfer function lowest critical frequency may be a pole. 
This situation is shown in Figure 4-8. One has, in effect, a realizable 
transfer function so long as all transfer function poles and zeros are 
restricted to lie along the negative real axis including the origin. A 
given transfer function can, therefore, be written as 

, N(s)/Q(s) , (4-26) 

Ej(s) D(s)/Q(s) 

where N(s)/Q(s) = Y (s) and D(s)/Q(s) = Y (s) and where Q(s) is 

A B 

selected so that both ^nd Yg(s) can be realized as two-terminal 

RC driving-point admittance networks. The synthesis of Y (s) and Y (s) 

A B 

may now be carried out individually. 

One method of synthesis is a partial-fraction expansion [Ref. 15] 

of Y (s) and Y (s) as indicated by Y(s) below: 

A B 
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Figure 4-8. Location of poles and zeros for 2-terminal RC network 
configuration 




Figure 4-9. Second Foster form realization of an RC admittance 
function 
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+ Kqo 



(4-2 7) 



Y(s) _ o , L 



+ 



s +<T 



s s 

Subsequent multiplication by s of each term in the right-hand side of 
equation (4-27) gives 



n 



Y(s) = 






K.s 

1 



s + 



+ K Qo s 



(4 -2 8a) 



i=l 



or 



n 






r . . 

Z_a -i_+ 



K 



00 



(4 -2 8b) 



i=l K. K.s 

1 1 



where K = 1/R , Kqd = Ccp and each l/K. = R. and each K./ (7~ . = C. 

O O ^ 11 l*^!! 

in units of ohms for the R's and farads for the C's. 

The resulting network for Y (s) and Y (s) is the configuration 

* A B 

shown in Figure 4-9. Note that this is the second Foster form for an RC 
admittance function . 

To give this synthesis method a more general application, Mathews 
and Seifert substituted three-terminal networks for the two-terminal 
networks, as shown in Figure 4-3. The transfer function of such a con- 
figuration was already developed in Section A of this chapter. The 
relation, assuming reciprocal networks , is 



E (s) 
o 

Ej(s) 






(4-29) 
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where and y (s) are transfer admittance functions. The trans- 

1 ZA i Z B 



fer admittance of a three-terminal RC network can have only simple poles 
[Ref. 16] which must lie on the negative real axis of the complex- 
frequency plane excluding the origin but may have zeros which lie any- 
where in the complex-frequency plane except on the positive real axis. 
Since these zeros need not be simple, the transfer function in equation 
(4-29) is perfectly general. 

The realization of a desired transfer function using the configuration 
in Figure 4-3 then boils down to synthesizing the transfer admittances in 
equation (4-29). 

Several procedures [Refs. 17-19] are available for synthesizing 
such transfer admittances but for an example, the method of Guillemin 
[Ref. 17] will be used. 

Example 4-2. Realize the high-pass transfer function derived in 
Example 2-3: 



The transfer function, H(s), in equation (4-29) may now be expressed 
as 



H(s) = - 



2 

s 



(4-30) 



+ 57.24 s + 5576.46 




(4-31) 



where 



yi 2 A^®) = - N(s)/Q(s) 



(4-32a) 



^12B^^^ = - D(s)/Q(s) 



(4 -3 2 b) 
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For y (s) and y (s) to be realizable, the roots of Q(s) must be 
12A 12B 

restricted to the negative real axis. Hence, Q(s) is judiciously ^ 
chosen to be 



Q(s) = (s + 50) (s + 150) 



such that 



y 




2 

s 

(s + 50) (s + 150) 



(4-33) 



(4 -3 4 a) 



, V _ (5576.46 + 57.24 s + s ) 

^12B (s + 50) (s + 150) 



(4 -3 4 b) 



By Guillemin's procedure, driving-point admittance functions are then 

chosen With the same poles as the transfer admittances, 

y (s) . These driving-point admittances are then expanded into RC 
1 2 B 

ladder networks in such a manner as to realize the zeros of the transfer 



admittances. For convenience, let 

, , . , _ (s + 10) (s + 97.5) 

^llA^®^ ^IIB^®^ (s + 50) (s + 150) 



(4-35) 



For reasons that will become apparent, y,«„(s) is then decomposed 

1 ZB 



as follows: 






(5576.46 + 57.24 s) 
(s + 50) (s + 150) 



(s + 50) (s + 150) ' 



(4-36) 



such that 



The choice is arbitrary but it can happen that for a particular 
choice, the spread in the calculated element values becomes too wide 
to be practical. 
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\ 



(4 -3 7a) 






(5576.46 + 57.24 s) 
(s + 50) (s + 150) 



''1262'=’ 



(s + 50) (s + 150) 



(4-3 7b) 



y ,„(s) is now expanded [Ref. 20] to realize the zeros of 

1 io IZd 1 

s = - 97.5 and at s = - CO : 



“'' 1161 '=’ 



28.5 



s + 97.5 



+ 1 + 



1 



(4-38) 



0.0156 s + 



6.4 



The resulting network, N , , is shown in Figure 4-10. 

B i 

Again y ,_(s) is expanded to realize the two zeros of V,ry^^(s) 
1 IB 1 Z dZ 

at s = 0: 






= 1 



7.7 



. (4-39) 



0.0109 s 



1.47 



1 



+ 6.3 



0.00158 s 



The resulting network, N , is shown in Figure 4-11. 

dZ 

Before these two networks, N and N , can be connected in 

B1 BZ 

parallel to realize y,ou(s) , it is first necessary to adjust their individual 

1 Z B 

admittance levels. From Figure 4-10, ~ ~ 1/7.693 = - 0.13 

at s = 0 and equation (4 -37a) gives yiou,(s) = - 0.746 at s = 0. A 

1 Z B 1 

multiplier for the network, is therefore 
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0.293 




Figure 4-10. 



Circuit realization of ^llBl^^^ 



0.0109 0.00158 6.3 




Figure 4-11 . 



Circuit realization of ^11B2^^^ 
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j 



I 

I 



0.174 



(4-40) 



K = 

B1 0.746 



For N , yioDO^®) “ “ 1/.63 = - 0.159 at s = - CD from Figure 4-11. 
dZ lZgZ 

From equation (4 -3 7b), “ “ 1 at s = - 00 . Consequently, 

1 Z BZ 






(4-41) 



Kg is then determined from the relation 



B 



K: 



1 1 



B1 



K 



B2 



(4-42) 



to be = 0.0832. 

B 

The admittance level of is now multiplied by 



K 



B 



K 



B1 



0.0832 

0.174 



= 0.478 



(4-43) 



and that of N is multiplied by 
BZ 



B2 



0.0832 

0.159 



= 0.523 



(4-44) 



The realization of observed to be exactly the network, 

N , except that the admittance level is multiplied by 
BZ 



K^ = Kb2 = 0.159 



(4-45) 



The complete network realization is shown in Figure 4-12. Note that 
the transfer function, H(s), in equation (4-30) is realized within the 
constant multiplier: 
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Figure 4-12 . Circuit realization of Example 4- 
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K 



A 



0.159 

0.0832 



K, 



1.91 



(4-4 6) 



B 



D. PAUL [1963] 

Paul [Ref. 21] sought to correct two limitations of the synthesis 
method of Mathews and Seifert: 

1 . No capability for independent adjustment of each transfer 
function coefficient, or, alternatively, of each pole and zero. 

2. Wide spread of element values. 

Using as a starting point the basic network configuration of Figure 

4-3 and its corresponding voltage transfer function, H(s) = - Yr,-,A^)/ 

y (s) , he proposed the adjustable network configuration shown in 
1 2B 

Figure 4-13 . 

It can be noted that the three basic T-networks, shown in Figure 
4-14 and governed by the relations in Table 4-2, are all integrated in 
both the input and the feedback networks . The configuration in Figure 
4-13 is a special case where 





'^2 = ' 




(4-47) 
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Figure 4-14. Type 1 T-Networks for Paul's configuration 
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Table 4-1. Some Transfer Functions Realizable by Figure 4-13 





V (s) 


Network Conditions 


' V^(s) 


All switches closed 


1 + k^sT + k^s^ T ^ 




1 + k sT +k^s^ r ^ 

4 6 






Switches 1 , 2 , 4 , 5 , 6 closed 


1 + k^s T 


Switch 3 open 


1 + k.sT + k-s^ T ^ 

4 6 


Switches 1,2, 4, 5 closed 


l+kjST 


Switches 3 , 6 open 


1 + k s T 
4 



Table 4-2 . Transfer Admittance Relationship of the 
Three Type 1 T -Networks 




where T = 2C^R^ = = 1/2C^R^ = CR 
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If 



C^=2C3=1/2C^ = C 
Rq = 1/2 



then 

T = CR . (4-49) 

Some transfer function forms possible for the configuration in 
Figure 4-13 are given in Table 4-1. 

Another special variation for the configuration shown in Figure 4-13 
is the case where 



Nj = l 



Nj = 2 



If 



”2 = = 2 






"C = '^B = '■'a = '" ' 



( 4 - 50 ) 



( 4 - 51 ) 



then 



= CR 



(4-52) 



In this case, the transfer function forms tabulatetd in Table 4-1 
result with every sT term in both numerator an(d (denominator replaced 
by 2sT. 
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If, instead of three standard T-networks, an RC configuration, such 
as the one shown in Figure 4-15, is used as both the input and feedback 
networks , the transfer admittance of each network can be shown to be 






(1 + sCRj) (1 + sCR^) 



[2 + sC(Rj + ^ ^4 



(4-53) 



+ [2 + sC (R 3 + R^)][Rj + 1^2 sCRjR 2 ] 



and the resulting voltage transfer function is 



H(s) 



(1 + sCR^) (1 + sCRj^) 
(1 + sCR„) (1 + sCR_) 

D U 



(1 + s T^) 0 + s Tp) 

(1 + S Tg) (1 + s Tg) ' 



where, in relation to the network in Figure 4-15, for the input network. 




R 



2 



= R 



B 



R 



3 




R 



4 




(4-55) 



and for the feedback network. 



I 





Figure 4-15. Type 2 T-Network for Paul's configuration 
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Figure 4-16. Type 3 Network for Paul's configuration 



87 







I 





i 



Note that independent adjustment of the time-constants, I , is possible 
by adjusting the relevant set of resistors and capacitors. 

If the network configuration, shown in Figure 4-16, is used as 
both the input and the feedback networks, the transfer admittance can 
be shown to be 



1 + sCR^ 



(4-57) 



and the overall active voltage transfer function is 



i + sCR^ 

1 + sCRg " ■ 1 +s Tg 



(4-58) 



In relation to Figure 4-16, for the input network, R = R and R = R ; 

i A ZB 



for the feedback network, R, = R„ and R„ = R. . 

IB 2 A 



The transfer function forms generated by the network configurations 
of Paul would ideally be used for realizing higher-order filter transfer 
functions, broken up into biquadratic polynomials and cascaded to form 
the desired transfer function. 

Two networks that would be especially useful for such realizations 
are shown in Figures 4-17 and 4-18. The network in Figure 4-17 is 
designed to simulate transfer functions of the form 



H(s) = - -- -Y ;p , (4-59) 

where T = RC . That in Figure 4-18 simulates transfer functions of the 
form 



88 
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Figure 4-17. Paul's circuit to simulate V (s)/V (s) = - 1/(1 + s T ) 

o 1 
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1 



H(s) = - 



1 



(4-60) 



1 + 2k. s 7* + s 7^ 

1 

where 1 . These two networks will be used in the example to follow. 

Example 4-3. Realize the fifth-order Butterworth transfer function 
of Example 2-2 . 

Reduced to a form suitable for simulation by Figures 4-17 and 4-18, 
the transfer function is 

H(s) = 1 . L_ 

1 + s 1 + 0.618 s + s^ 

62.83 62.83 (62.83)^ 

. 1 . (4-61) 



1 + 1.618 s + 
62.83 



2 

s 



(62.83)^ 



Comparison of the first product term with equation (4-59) yields 
the relation 



T = 1/62.83 



(4-62) 



For C = 1 , the corresponding value of R is determined to be 

.6 



R = 



10 



= 15.9 K il 



(4-63) 



62.83 

Comparing the second and third product terms with equation (4-60) 



gives the following potentiometer settings: 
kj = 0.618/2 = 0.309 
k^ = 1.618/2 = 0.809 



(4-64) 

(4-65) 



The final circuit realization is shown in Figure 4-19. 
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Figure 4-19. Circuit realization of Example 4- 



E. HAKIM [1964] 

Hakim [Ref. 22] considered the realization of minimum phase 
transfer functions. His method essentially follows the same line of 
approach as that of Mathews and Seifert. His method of development, 
however, is more detailed and gives a clearer picture of the synthesis 
procedure . 

If the general transfer function is 

H(s) = - N(s)/D(s) 

_L _L 

as +a ,s +....+a,s+a 
_ _jn m-1 1 o_ 

, n , n-1 

bs +b ,s + ....+b,s+b 

n n-1 1 o 



m 



(s + c>^,) (s + ^ ) . . . (s + cxC ) 
. 1 z m 



^n ^ ^ . . . . (s + 



(4-66) 



where a to a and b to b are all positive real coefficients and where 
o m on 

the negative sign has been introduced in anticipation of the phase in- 
version caused by operational amplifier action, for the transfer function 
to be physically realizable, the order m of the numerator polynomial N(s) 
must be less than or equal to the order n of the denominator polynomial 
D(s) . Stability considerations restrict all the poles ~ 
left half of the complex-frequency plane. If all the zeros - to 
~ C>C^ also restricted to the left half of the plane, then the transfer 
function is said to be minimum phase. 
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For the feedback network configuration in Figure 4-3, the open- 



circuit voltage transfer function has been shown to be 



H(s) = 



E (s) 
o 

Ej(s) 






(4-67) 



If the transfer admittance functions, and y,oD(s) are to be 

z lA 1 ZB 

realized by RC networks, an auxiliary polynomial Q(s) has to be chosen 

with roots restricted to the negative real axis such that Q(s) = 

(s + ^ ) (s + .... {s + ) where ^ to ^ are real positive 

j. z o j. n 

numbers . 

Equation (4-66) may be written as 



H(s) = 



N (s)Z.Q.(s). 

D(s)/Q(s) 



and, equating equations (4-67) and (4-68), 



(4-68) 



y2iA(s) = - N(s)/Q(s) 

= - D(s)/Q(s) . (4-69) 

Assuming the transfer function, H(s), has a pair of complex conjugate 

poles or zeros, y (s) and can be realized using Guillemin’s 

21A 12B 

procedure [Ref. 17]. Equation (4-69) may now be written as 



rn ^ ni-1 ^ ^ ^ 

as +a ,s +....+a,s+a 
m m-1 1 o 



^21A^®^ ■ (s + cT^) (s + 



(s t ^^) 



(4-70a) 
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+ .... + b, s + b 
1 o 



(4-70b) 



, n , , n-1 

b s + b , s 

( \ _ _n Ilzi 

''i2b'^ ' (s + cTy) (s + 



(s + cr^) 



Following Guillemin's procedure, short-circuit driving-point admittances 
are arbitrarily chosen. For convenience, let 



(s + (^ j) (s + (^ 2 ) • ♦ ♦ (s + ($ j^) 

(s) = y, ,„(s) = ; — . ■■>» , — . ^ ■ v , (4-7 

22A IIB (s + ^) (s + ... (s + ) 

1 Z n 

where , to X are all positive numbers. Since N. and N„ are to be 
^ 1 ^ n A B 

RC networks, the poles and zeros of equation (4-71) must alternate and 
the lowest critical frequency must be a zero. 

Synthesis of each network may now follow. For y21A^^^' 
first decomposed such that 



N(s) = n^(s) + n 2 (s) = . . . + n^(s) , (4-72a) 

where 

n, (s) = a + a, s 
1 o 1 

f \ 2 ^ 3 

tajS 



n (s) = a ^ + a s*' . (4-72b) 

r n-1 n 

The last polynomial n^(s) contains one or two terms depending on whether 
the order n of N(s) is odd or even. 
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V21a‘=>= 



A '' 2 IAI*®’ ''A ''21A2*®* ■'' 



K 



A2 



•^A ''21Ar‘=> 



K 



Ar 



(4-75) 



and if the resulting scaled ladder realizations of each of these components 
are connected in parallel at their input and output ports, the resulting 
two-port network will realize the short-circuit admittance 7^22^^^^ 
the required zeros of • 

The transfer admittance y,-,r,(s) is realized in a similar manner. 

The resulting network is shown in Figure 4-20. 

Although basically the same method has been used in Section C to 
realize a high-pass transfer function, another example using Hakim's 
method to realize a low-pass filter will be given for purposes of com- 
parison with other low-pass realizations. 

Example 4-4. Realize the Butterworth low-pass transfer function 
of Example 2-1 . 

The transfer function 



H(s) = - 



3947.84 



s + 88.84 s + 3947.84 



(4-76) 



is first divided in both numerator and denominator by a judiciously 
chosen Q(s) = (s + 30) (s + 60) so that 



( \ = 3947.84 

■ (s + 30) (s + 60) 



(4-77a) 
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and 



3947.84 + 88.84 s 
(s + 30) (s + 60) 



(4 -7 7b) 



^128 



(s) = - 



2 

s 

(s + 30) (s + 60) ' 



where 



. ^ 3947.84 + 99.84 s 

^12B1 ® ” (s + 30) (s + 60) 



(4-78a) 



and 



^12B2 



(s) = - 



2 

s 

(s + 30) (s + 60) 



(4 -7 8b) 



Driving-point admittances are then chosen. For convenience, let 



V (s) = V (s) = ♦ (4-791 

^22a' ' ^IIB^ ^ (s + 30) (s + 60) ’ 

(s) is first expanded into a ladder network in such a manner as to 

(s) at s = - 44.5 and at s = - 00 : 



"IIB 

realize the zeros of y 



12B1 






1 

6.51 + 1 + 1 

s +44.5 s + 1 

28.99 2.899 



(4-80) 



The network realization, N . , is shown in Figure 4-21. 

Bi 

y ,_(s) is next expanded to realize the two zeros of y,oDo^^^ 
1 1 B iZBZ 

at s = 0; 
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o 



0.1465 

rVVV^ 



1 2.899 

~V\AA-t — v^AA 



0.1536 



0.0345 



N 



B1 



Figure 4—21. Circuit realization of end Y ^ 




Figure 4 22 . Circuj*; realization of y (s) and y (s) 

iZBZ 1 1B2 
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1 



4.04 



•(4-81) 



0.01795s 



+ 



1 .49 



+ 12.35 



1.686s 



The network realization, N „ , is shown in Figure 4-22 

BZ 



Scaling factors are determined to be 



^B ^ 0.0471 

0.1127 
B 1 

and 



0.418 



(4-82a) 



B _ 0.0471 



B2 



0.081 



= 0.582 



(4-82b) 



y (s) is then expanded to realize the double zeros of y (s) at 
s = - 00 : 



^22A<=> = 



(4-83) 



1 + 



0.0282s + 



2.18 + 



0.0423s + 



0.865 



The scaling factor for the input network, N , is calculated to be 



^ 0.247 
A 2.195 



0.1125 



(4-84) 



The final circuit realization is shown in Figure 4-23. 

The transfer function in equation (4-76) is realized within a 
constant multiplier 



^ = 0.1125 = 2.39 

K_ 0.0471 

" 100 



(4-85) 



I 

I. 



I 
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Figure 4-23. Circuit realization of Example 4- 



F. MUIR AND ROBINSON [1968] 

One method proposed by Muir and Robinson [Ref. 23] utilizes 
three operational amplifiers to realize a second-order transfer function. 

A three-amplifier realization of a low-pass filter is shown in Figure 4-24. 
The basic building block for this configuration is a single operational 
amplifier feedback circuit, shown in Figure 4-25, with the familair 
voltage transfer function, H(s) = - Y^/Y^. 

If the currents at the input of each operational amplifier of 
Figure 4-24 are added, assuming A , the resulting equations are 

''l ^6 '"o ^3 ^ ^4 = ^ 

\ ^1 ^^2 ^^2^ " ® (4-86b) 

g + V sC = 0 . (4-86c) 

B 5 o 1 

Solving the above equations simultaneously yields the following 
voltage transfer function: 



V 






s + (g/c^) S + (SiSjS/g^CiCj) 



(4-87) 



If equation (4-87) is identified with the prototype low-pass filter trans- 
fer function in equation (2-16) , 



V 

o 




H u> 



2 2 
s + ^ tc> s + uJ 
o o 



(4-88) 
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Figure 4-25. Basic circuit block used in triple ampli'ier realization 
of transfer functions 
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B 





it can be seen that 



H ^ 

0 




(4-89a) 


a 




(4 -89b) 


2 

^0 


= 3 , 5395 / 940^03 . 


(4-89c) 



To realize a desired low-pass transfer function, it is only necessary to 
make the following steps: 

1. Select Cj, C 2 and g^. 

2. Determine g^ from equation (4-89b). 

3 . Let 93 = • 

4. Employ the relation ^1 _ ^5 , . , 

— - — = CO^ to obtain g^ and g^ 

2 1 

5. Finally compute g from equation (4-89a). 

b 

High-pass, bandpass and band-reject configurations, similar to Figure 
4-24 may be obtained from Ref. 23. 

Note that in the realization of higher-order transfer functions, 
the configuration of Figure 4-24 reverts to the status of a second-order 
building block . 

Example 4-5. Realize the low-pass Butterworth filter transfer 
function of Example 2-1. 

Following the design procedure steps laid out in this section, 
let 

Cj = = 1 yt^F (4-90a) 

R^ = R = 100 kH . (4-90b) 

3 4 
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From the relation, ~ ^8.84, is calculated to be 



= 11.27 K£1 



(4-91) 



Knowing CO = 62.83 rad/sec, R. and R are determined to be 

O X ^ 




(4-92) 



Finally, R_ is calculated from equation (4-89a) to be 

D 





H CO 



1 



- = R^ = 100 kH 



(4-93) 



o 



since H = 1 . 

The final circuit realization is shown in Figure 4-26. 

G. ADVANTAGES AND DISADVANTAGES 

Compared to comparable realizations falling under the other 
categories in this paper, single-feedback realizations will be observed 
to require considerably more elements . This can be noted especially 
in the method of Hakim, where the number of ladder networks increases 
with the order of the transfer function. This disadvantage is due to 
having to realize the input and the feedback networks separately, each 
network with its own set of natural frequencies. Because of this same 
reason, however, such realizations determine transfer function pole 
locations by the zeros of the two passive networks, rather than by the 
gain of the operational amplifier, and hence, are quite stable. 
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Figure 4-26. Circuit realization of Example 4- 



V. INFINITE-GAIN MULTIPLE -FEEDBACK REALIZATIONS 



This chapter contains the significant contributions falling under 
the general category of Infinite-Gain Multiple-Feedback Realizations. 

The difference between this category and that in Chapter IV is that, 
as the classification implies, the feedback connections are made to 
more than one point in the input network from the output of the operational 
amplifier. Realizations falling under this category will thus be of the 
general network configuration shown in Figure 5-1. 

A. THEORETICAL DEVELOPMENT 

Consider the network configuration shown in Figure 5-2. This is 
the second -order prototype of the general multiple-feedback configuration 
in Figure 5-1. Consisting of five two-terminal admittances and one 
infinite-gain operational amplifier, this network would be an ideal choice 
for illustrating Nathan's method [Ref. 12] for analyzing networks con- 
strained by an operational amplifier. 

The Y matrix for the network without the operational amplifier is 



1 


2 


3 


5 








0 


0 


1 


-^1 






-''4 


2 

(5-1) 


0 


-^3 




-^5 


3 


0 


-^4 




(Y 4 +Y 5 ) 


4 
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Figure 5-1. General infinite-gain multiple-feedback configuration 
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Figure 5-2. Second-order prototype of multiple-feedback configuration 
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In accordance with Nathan's procedure, the desired matrix, [Y‘] , for 
the network with the operational amplifier is obtained by deleting the 
fourth row and the third column from the Y matrix in equation (5-1): 



1 2 4 





-Y. 



(Y +Y +Y +Y ) 
' 1 2 3 a ' 



-Y. 




(5-2) 



To obtain the voltage transfer function, E^(s)/E^(s), equation (C-7) in 
Appendix C is applied: 



£4(3) y ^4 

E^(s) 



-Y, 






-Y, 



(Y «+Y +Y +Y ) 
^ 1 2 3 a ’ 



-Y. 



-Y. 



-Y, 



(5-3a) 



or 



H(s) = 



-Y Y 
13 



Y (Y +Y +Y +Y ) + Y Y 
5 1 2 3 a ' 3 4 



(5-3b) 



B. RAUCH AND NICHOLS [1956] 



Using equation (5-3) as the basic equation, Rauch and Nichols 
[Ref. 24] developed what was to be referred to as the "Rauch" filter. 
The filter network, shown in Figure 5-3, is the low-pass version of the 
general second-order prototype in Figure 5-2. 



no 



The transfer function corresponding to the network in Figure 5-3 



is 



H(s) = - 



2 

sCC +sC (G +G +G) + GG 
2 5 5 1 3 4 3 4 



(5-4a) 



or, in a more convenient form 



1 



H (s) = - 






(5 -4b) 





1 



The synthesis procedure then becomes a matter of equating the 
coefficients of the desired second-order low-pass transfer function 
of the form of equation (2-16) to the element coefficients of equation 
(5-4b). 

An example of this procedure is given in Example 5-1 . 

Example 5-1. Realize the low-pass transfer function in Example 
2-1 using a second-order Rauch filter. 

The transfer function, repeated here for convenience, is 



H(s) = - 



3947.84 



(5-5) 



s^ + 88.84 s + 3947.84 



Equating coefficients of equations (5-4b) and (5-5) yields 




88.84 



(5-6a) 



and 



I/C 2 C 5 R 3 R 4 = 3947.84 
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(5-6b) 





Fj .e 5-5. High-pass active network configuration of Bridgman 
anw Brennan 
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H(s) = 




(5-10) 




are capacitors as shown in Figure 5-6a, the transfer function is 



Again element values for a desired second-order transfer function are 
determined using the coefficient-matching approach. 

Example 5-2. Realize the Chebyshev high-pass transfer function 
of Example 2-3 . 

The transfer function is converted to the form 



H(s) = 




(5-11) 



If Z and Z are capacitors as shown in Figure 5-6b, the corresponding 
transfer function is 



H(s) = 




. ( 5 - 12 ) 



H(s) = 




(5-13) 



2 

(s /5576.46) + (57.24/5576.46) s + 1 



Equating coefficients of equations (5-10) and (5-13) gives 



C C R R = 1/5576.46 

J. o ^ O 



(5-14a) 



C C R R = 1/5576.46 

O T ^ O 



(5-14b) 



R (C + C + C ) = 57.24/5576.46 

Zt 1. o ^ 



(5-14c) 
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Figure 5-6b. Bandpass configuration No. 2 of Bridgman and Brennan 
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H(s) = - 



(5-18) 



a-s +a,s + a 
2 , i o 

b + b + b s + 1 

O ^ i. 

where all numerator and denominator polynomial coefficients are positive 
and real . 

Wadhwa investigated the realization of three separate cases of 
equation (5-18). The first of these [Ref. 26] is the realization of the 
case where a^ = =0 such that 

a 

H(s) = - ^ 1 • (5-19) 

b^s + b^s + bj^s + 1 

Three possible arrangements of resistors and capacitors to realize 
equation (5-19) are shown in Ref. 2 6 with their corresponding relations 
between transfer function coefficients and network elements. One such 
arrangement with its corresponding coefficient -component relations is 
shown in Figure 5-8. This particular arrangement is an extension of the 
second-order Rauch filter. 

For the realization of the special case where 



H(s) 





3 2 

b s + b s + b s + 1 

O u <i- 



(5-20) 



three possible combinations of resistors and capacitors are shown in 
Ref. 27. 

One combination which is an extension of the network in Figure 
5-6a is shown in Figure 5-9. Its coefficient -component relations have 
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i 





Mr 




^0 = 



bj = 1/3 [(5bg+l)RCg + (bjj+l)RCg] 
b^ = 1/3 RCg Ubj^RC^ + (2bj|+l)RCg] 



^=>’0 






^7 = 



b s. b (b +1) 
1 > ^ 0 

b„ 



Figure 5-8. Wadhwa's configuration to simulate 
H(s) = - aQ/(b 2 S^ + b^s^ + b^s + 1, 
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Figure 5-9. Wadhwa's active network configuration to simulate 
H(s) = - a^s/Cb^s^ + + b^s + 1) 




Figure 5-10. V/adhwa's active network configuration to simulate 
H(s) = - a^s^/ih^s^ + b^s^ + b^^s + 1) 
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been intentionally omitted because the realization conditions become 
too complicated to be of any illustrative value. 

Finally, Ref. 2 8 treats the realization of the special case where 



H(s) = - 



^2" 

3 2 

b^s +bj^s + l 



(5-21) 



One of the three possible combinations discussed is shown in Figure 
5-10. 



Realizations of still higher-order transfer functions generalized 



by 



H(s) = - 



, n-1 ^ n-2 . . , 

a (a ,s +a „s +...+a,s+l) 
o n-1 n-2 1 

b s^ + b ,s^ ^ + ... +b_s^+b,s + l 
n n-1 2 1 



(5-22) 



following the same approach of coefficient-matching subject to certain 
realization conditions have been described in Ref. 29. The details of 



these, however, will not be included in this paper. 



E. FOSTER [1965] 

Foster [Ref. 30] developed a method of realizing low-pass filter 
transfer functions of orders higher than the third by using Rauch's two- 
and three-pole filter configurations . 

Basically, if a fourth -order low-pass filter is desired, the cor- 
responding transfer function is realized by cascading two two-pole 
Rauch filters. If a fifth -order filter is desired, a three-pole and a two- 
pole filters are cascaded. 
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In general, if an odd -order low-pass filter is desired, a three- 
pole Rauch filter followed by as many two-pole filters as necessary 
effects the realization. If an even-order filter is required, the necessary 
number of two-pole Rauch filters is used. 

The two-pole or the second -order Rauch filter has already been 
developed in Section B of this chapter. The third -order Rauch filter is 
shown in Figure 5-11. Its transfer function is 



E (s) 
o 

E^(s) 






(5-23) 






'^ 3 '' 5^'=4 















’^l'^3'^5'*7^2^4''6 



A third-order low-pass filter may then be realized by equating the 
coefficients of its transfer function to the coefficients of equation (5-17). 
It is normally convenient in this kind of procedure to make all resistors 
of the same value, choosing the resistor value to come up with practical 



capacitor values. 

Using equal resistors for both two- and three-pole Rauch config- 



urations, Foster simplified the realization of higher-order low-pass 
filters by tabulating normalized capacitor values for Butterworth, Bessel, 



and Chebyshev filters up to order 10. A portion of the capacitor table 
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"■ f 

Figure 5-11. Third-order Rauch active filter 
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for Butterworth filters is shown in Table 5-1. The values in the tables 



are expressed in farads. To denormalize, each capacitor value is 

divided by R ^ , where R is the chosen resistor value and U) is the 
o 

angular cutoff frequency. A good practice in choosing R is to check if 

dividing the largest and the smallest capacitor values by R results 

o 

in practical values. 

To illustrate this method, an example follows. 

Example 5-3. Realize the fifth-order Butterworth filter transfer 
function in Example 2-2. 

The angular cutoff frequency of the desired filter is = 62.83 
rad/sec. From Table 5-1, it is noted that the largest capacitor value 
is 4.31481 farads and the smallest is 0.21386 farad. An appropriate 
value for R is therefore 10 kH- so that ^ = 6.283 x 10^. The required 

capacitor values are therefore: 



2.16471 

- c oar, ^.45 UF 

1 6.283 X 10 / 



= e.87„F 

6.283 X 10 






0.21386 
6.283 X lo' 



0.340 yUT 



^4 = 



1.85410 
6.283 X lo' 



= 2.960 



(5-24) 
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Table 5-1. Normalized Capacitor Values for Butterworth Filter 



n 


<=i 


C 

2 


=3 




^5 


1 


1.00000 










. 2 


2.12132 


0.47140 








3 


2.37484 


2.59100 


0.32503 






4 


3.91969 


0.25512 


1.62359 


0.61592 




5 


2.16741 


4.31481 


0.21386 


1.85410 


0.53935 
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0.53935 



0.859 JaY 



C 



5 



6.283 X 10 



5 



The resulting filter network is shown in Figure 5-12 . 



F. HOLT AND SEWELL [1965] 

The method of design proposed by Holt and Sewell [Ref. 31] 
differs somewhat from the basic multiple-feedback configuration shown 
in Figure 5-1 in that it involves two cascaded sections shown in Figure 
5-13. However, the active section includes an infinite-gain operational 
amplifier and^a multiloop feedback circuit, as shown in Figure 5-14, and 
hence, belongs under the multiple-feedback classification. 

The method is limited to the realization of biquadratic functions 
having transmission zeros on the imaginary axis of the s-plane. Such 
functions are known as Cauer or elliptic functions, having low-pass 
transfer functions of the form: 



2 

s + a 



H(s) = 



s + b, s + b 
1 o 



(5-25) 



The active circuit of Figure 5-14 is used to produce the transfer function 



s + A,s + A 

H (s)= -2 ^ ^ 

^ s + b,s + b 

1 o 



(5-26) 



and the passive section is used to produce the transfer function 
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Figure 5-14, Active section of Holt and Sewell 




Figure 5-15. Active network of Holt and Sewell to realize 

/ X 2 

H (s) = s + A^.s + Aq 
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H2(s) = 



(5-27) 



2 

s + a 
o 

+ A, s + A 
1 o 

The active and passive sections are then cascaded to yield the required 
transfer function in equation (5-2 5). 

Analysis of the operational amplifier section yields the transfer 
function 



YY +Y (Y +Y +Y) 

H fs) = _ -i-2 Iz V 

1 Y^Y^ +Y, (Y, +Y^ +Y^) 



(5-28) 



2 3 



1 



One possible combination of elements that would give the desired form 
of equation (5-26) is for Y^ and Y^ to be resistors, Y^ and Y^ capacitors 
and Y to be a resistor and a capacitor in parallel, as shown in Figure 

D 

5-15. The transfer function in terms of the elements is 






R R 

H(s) = - ^ 



R, 



R. 






R1R4 



R. 



(5-29) 



Equating coefficients of equations (5-26) and (5-29) gives: 



b = 1/R.R. 
o 14 



h ' '=2 " 



1 



(5-30) 



If 






(5-31) 
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the resulting resistor values are 



R, = b ,/2b 
1 1 o 

= 2/b^ . (5-32) 

If R = 1 , equation (5-29) becomes 
o 



Hj(s) = 



2b 



+ 1 \ s + s 



(5-33) 



b + b, s + s 
o 1 



The passive network following the active section has to produce the 
following transfer function: 



s^ +a 



H2(s) = 



2 + / 2b 



s + b 



(5-34) 



Equation (5-34) may be synthesized following Guillemin's parallel- 
ladder method [Ref. 20] provided it has only real poles, or expressed 
in terms of its denominator coefficients, provided that 




(5-35) 



If the desired transfer function of the form of equation (5-25) satisfies 
this realization condition, the cascaded sections of Holt and Sewell 
may be used . 
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Example 5-4. Realize the first product term of the transfer function 
in Example 2-6 . 

The product term 



. s^ + 3908.36 (5-36) 

ms; — 2 

s + 5.51 s + 3420.98 

is first checked to determine if it satisfies the realization condition of 

equation (5-35) . With b = 3420 .98 and b, =5.51, it can be seen that 

o 1 

equation (5-3 6) does satisfy the realization condition. 

From equation (5-33), the active section transfer function is 



H^(s) = 



2( s + 1243 s + 621) 
s^ + 5.51s + 3420.98 



(5-3 7) 



The passive section transfer function, from equation (5-34), is 



2 

„ / , s + 3908.35 , (5-38) 

n is; - - 2 

s + 1243 s + 621 

where a = 3908.36. 

0 

Equations (5-33) and (5-34) were determined by setting C_ = C„ 

= C = 1 and R = 1 in the active section, so that all that remains to 

O D 

specify the active section completely is to determine R^ and Rj^ as 
follows: 

R=2/b= 0.363 (5-39a) 

4 1 

R = by2b = 0.805 x 10"^ , (5-39b) 

1 i o 
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where C's are in farads and R's in ohms. 



The synthesis of the passive transfer function, H 2 (s), involves 
choosing an appropriate Q(s) . In this case, one choice is 

Q(s) = (s + 2) (s + 1500) , (5-40) 

such that 



s^ + 3908.36 

(s + 2) (s + 1500) 

s + 1243 s + 62 1 
(s +2) (s + 1500) 



(5-41) 



where 

2 

( \ = s + 3908.36 (5-42a) 

^12P'^^ “ (s + 2) (s + 1500) 

. . ^ s^ + 1243 s + 621 . (5-42b) 

^22P'^^ (s + 2) (s + 1500) 

y (s) can now be expanded so as to realize the double zeros of 

^ L 

yi 2 pi(s) = - sV(s +2) (s + 1500) at s = 0: 



^22Pl("^ = 



1 

4.83 



1 

0.310s 



1.61 



1 

0.000115s 



+ 5.81 



y 22 p(s) can then be expanded to realize the double zeros of 
yi2p2(s) = - 3908.36/(s + 2 ) (s + 1500) at s = - CO : 
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Y 

22P2 



(s) = 



1 + 



1 

1 

s + 1 

259 0.21 + 1_ 

0.55s + 



1 ^ 

3.62 



(5-44) 



To be able to connect the two ladder-network components of y (s) , 

X ^ 1 

shown in Figure 5-16, their admittance levels have to be individually 
adjusted by the following scaling factors: 



S ^ 0.0276 

Kp^ 0.207 



0.1333 



(5-45a) 



^ ^ 0.0276 

Kp2 0.0329 



0.839 



(5 -4 5b) 



The complete normalized circuit realization is shown in Figure 5-17. 
Both the active and the passive sections may be independently de- 
normalized to give practical element values. 



G. ADVANTAGES AND DISADVANTAGES 

For second-order functions, this category of realizations has the 
distinct advantage of usually requiring a minimum number of elements. 

For instance, in the realization of the low-pass transfer function of 
Example 2-1, Hakim's method of realization, belonging to the single- 
feedback category, required one operational amplifier, nine resistors 
and six capacitors. The Rauch filter, belonging to this category, re- 
quired one operational amplifier, three resistors and two capacitors. 

For higher-order cases, Wadhwa and Aggarwal have shown that realization 
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0.000115 0.310 




Figure 5-16a. Network realization of ^22P1^^^ 




Figure 5-16b. Network realization of ^22P2^^^ 
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Figure 5-17. Circuit realization of Example 
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i 



of such functions is possible by simply extending the basic second-order 



prototype configuration but the determination of the element values 
becomes increasingly complicated with increasing order because of the 
nonlinear nature of the relating equations . Such higher-order cases can 
easily be realized by factoring them into first- and second-order product 
terms and cascading the realizations of these terms. Cascading terms 
is possible without appreciable interaction because the output impedance 
of each realization is that of the operational amplifier. 
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VI. INFINITE -GAIN GENERAL REALIZATION TECHNIQUES 



Two general realization procedures using infinite-gain operational 
amplifiers are presented in this chapter. They are general in the sense 
that they can realize any rational function of any order. As opposed to 
the realization techniques of the last two chapters, these two procedures 
are not limited by the restrictions on the passive RC configurations that 
make up their input and feedback networks. 

A. LOVERING [1965] 

The general circuit proposed by Lovering [Ref. 32] is shown in 
Figure 6-1. Since the two operational amplifiers are assumed to have 
gain approaching infinity, the current summation at each of their summing 
points may be written directly as: 




(6-1) 



V Y„ + V Y, + V Y,, = 0 
12 5 o 6 



(6-2) 



If Y =Yj-, subtracting equation (6-2) from equation (6-1) yields 

o O 




(6-3) 



Putting equation (6-3) in transfer function form gives 



V 




o 



H(s) 



(6-4) 



V 



1 



Since the general transfer function may be expressed in the form 
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H(s) = Nls} 
D(s) 



(6-5) 



if both numerator and denominator are divided by a polynomial, 0(s), 
having simple roots restricted to the negative real axis and whose degree 
is one less than that of N(s) or D(s), whichever is greater, the transfer 
function may then be written as 



H(s) 



N(s)/Q(s) 

D(s)/Q(s) 



If N(s)/0(s) is expanded [Ref. 15] in the partial-fraction form 



( 6 - 6 ) 



n 



N( s ) ^ j. 
D(s) o 



i=l 



K.s 

^OOs 



or, expressed in terms of R's and C's, 



n 



(6-7a) 



= r * ) _J_ \ ^ ^ Ceos . (6-7b) 

° R. C.s 

i=l 1 1 

the result is a string of positive and negative terms because the 
vatious K's are all real but may be positive or negative. If 
in equation (6-4) is now equated to N(s)/Q(s), Y^ can directly be 
associated with all the positive terms in equation (6-7) and Y 2 with all 
the negative terms. 

A similar expansion of D(s)/Q(s) can then be carried out and all 
the positive terms of the resulting expansion can be equated to Yg and 
all the negative terms to Y^. 
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An example of this procedure will clarify the steps indicated . 
Example 6-1. Realize the low-pass Butterworth transfer function 
of Example 2-1: 



H(s) = 



3947.84 



s + 88.84 s + 3947.84 



=' N(s) 
D(s) 



(6-8) 



If the dividing polynomial is chosen to be Q(s) = s + 100, the resulting 
numerator rational polynomial can be expanded as 



H(s) ^ .Q 39.48s . (6-9) 

Q(s) " s + 100 

Similarly, the denominator rational polynomial can be expanded to 
produce 

• ^ = 39.48 

Q(s) s + 100 



Equating equation (6-9) to Y 2 and equation (6-10) to Yg - the 
following relations can be made: 



Y^ = 39.48 
39.48s 

^2 s + 100 
Yg = 39.48 + s 

50.64s 

M s + 100 



(6-1 la) 
(6-llb) 
(6-llc) 
(6-lld) 



The following element values can then be calculated: 
Rj = 0.0253 ohm 

R„ = 0.0253 ohm 
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0.3948 farad 



C 



2 



R 



6 



0.0253 ohm 



C 



6 



1 farad 



R 



4 



0.01971 ohm 



C ^ = 0 . 5064 farad 
4 



(6-12) 



5 



If a scaling factor of 10 is applied to the values in equation (6-12), 
the following component values result: 

R^ = 2.53 Kil 

R^ = 2.53 KH 
= 3.948 yWF 

R^ = 2.53 Kjn 
b 

Cj. = 10/.F 

R, = 1.971 Kil 
4 

= 5.064 y«F . (6-13) 

R* and R may conveniently be chosen to be 1 K IT . 

O 0 

The complete realization is shown in Figure 6-2 . 

To illustrate the property that this procedure can indeed realize 
any transfer function, even with poles and zeros in the right half of 
the complex-frequency plane, consider the next example. 

Example 6-2 . Realize the transfer function 
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H(s) = 



2 

s - s + 4 



(6-14) 



NM 



2 

s + 12s + 2 



D(s) 



If Q(s) = s + 2 is chosen 



Nis) 




(6-15a) 



Q(s) 



and 







(6-15b) 



Q(s) 



If and are chosen to be 1 ohm, the normalized circuit realization 



is shown in Figure 6-3. 

B. BRUGLER [1966] 

Brugler [Ref. 33] used a differential-input operational amplifier 
to develop his general realization technique. The configuration he used 
is shown in Figure 6-4. 

A current summation at the two input nodes of the operational 
amplifier gives: 




(6-16) 



and 




(6-17) 



Solving equations (6-16) and (6-17) in terms of V and V results in: 

A B 



V 



A 




(6-18) 
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Figure 6-4. Brugler's active network configuration 



and 



''b' 






(6-19) 



But since an infinite— gain operational amplifier forces its two inputs to 
be equal, 

Va=V . (6-2 0) 



So, equating equations (6-18) and (6-19) gives the transfer function of 
the circuit in Figure 6-4: 



V 



o 




^1 '^4 ^ y.S - ^4 <^1 ^ ^2 

’'e <’'i + ^2 V - ""a '''4^^5-^V 



If a simplifying condition is now applied. 



’'i ^ ""a ’^3 = '^4 * ^5 



the transfer function in equation ( 6 - 9 ) becomes 



( 6 - 21 ) 



( 6 - 22 ) 



V 



o 





(6-23) 



Equation (6-23) is exactly the form of the transfer function of 

Lovering's general circuit and exactly the same realization techniques 

can be applied to generate the required elements. 

With Y^, , Yg , and Y^ known, equation (6-22) can now be used 

to determine Y and Y by setting 
^ o 

^2 = Y 4 + Y 6 , (6-24a) 

Y = Y + Y 
5 1 ^3 
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(6-24b) 



Common terms in Y, + and Y, + Y^ can be subtracted out without 

4 6 13 

violating equation (6-22), thereby simplifying the resulting realizations 
of Y 2 and Y^ . 

An example will illustrate this procedure. 

Example 6-3. Realize the Butterworth low-pass transfer function 
of Example 2-1 . 

This same transfer function has already been realized using 
Lovering’s method in the previous section and if the same partial- 
fraction expansion is used for a chosen Q(s) = s + 100, the following 
relations can be made: 



39.48 


(6-25a) 


39.48s 
s + 100 


(6-25b) 


39.48 + s 


(6-25c) 


50.64s 
s + 100 


(6-25d) 



If equations (6-24) are now applied, Y 2 and Y^ are determined to be 



Y 



2 



39.48s 
s + 100 



+ 39. 48 + s 



(6-26a) 



Y 



5 



39.48 + 



50.64s 
s + 100 



(6-26b) 



Note that 39.48 is a common term in equations (6-2 6a) and (6-2 6b) and 
may be subtracted out. 

If a scaling factor of 10^ is now used, the final realization of 
the desired transfer function is shown in Figure 6-5. 
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2.53 Kil 
-AA^ 



I 



1.971 KiZ 5.064 M F 



2.53KH' 



10 /X F 



3.948 
/< F 



J 



2.53KG 3.948Ar_ 
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f 1.971 Kil 
!:z 5.064 /KF 




10 /<? 



Figure 6-5. Circuit realization of Example 6-3 




Figure 6-6. Circuit realization of Example 6-4 
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I 






Example 6-4. Realize the bandpass transfer function of Example 



2-4: 



H(s) = 



62.83 s 

s^ + 1.257 s + 3947.84 



msi 

D(s) 



If Q(s) = s + 100 is chosen, the expanded terms are 



(6-2 7) 



NM 

Q(s) 

and 

D(s) 

Q(s) 



62.83s 
s + 100 



39.48 + s 



By association. 



62.83s 
s + 100 

0 

39.48 + s 

138.22s 
s + 100 



138.22s 
s + 100 



(6-28a) 
(6-2 8b) 

(6-29a) 

(6-29b) 

(6-29c) 

(6-29d) 



Use of equation (6-24) gives 



= 39.48 + s 

62.83s + 138.22s 

5 s + 100 s + 100 

From equations (6-29) and (6-30), 

determined: 



(6-30a) 

= 211.05s . (6-30b) 

s + 100 

the following element values are 
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1/62.83 = 0.0159 



Cj = 62.83/100 = 0.6283 
= 1/39.48 = 0.0253 

C2= 1 

R = 1/138.22 = 0.00723 

C = 138.22/100 = 1.3822 

R = 1/211.05 = 0.00474 
o 

C = 211.05/100 = 2.1105 
o 

R^ = 1/39.48 = 0.0253 
6 

= 6 = * 



( 6 - 31 ) 



where R's are in ohms an(d C's in faracds. 

0 

If a scaling factor of 10 is applieid to equation (6-31), the 
resulting circuit realization is shown in Figure 6-6. 



C . ADVANTAGES AND DISADVANTAGES 

The two procedures discussed in this chapter retain all the 
advantages associated with techniques using infinite-gain operational 
amplifiers. The realized circuits are relatively less sensitive to the 
variations in gain of the active element. Cascading of low-order 
realizations can be accomplished without appreciable interaction between 
them because of the operational amplifier's very low output impedance. 
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And, in addition, these two procedures are capable of realizing any 
transfer function that a designer might feel is desired. It must be 
appreciated, however, that the realizations are not necessarily optimum 
as far as number of required elements is concerned. 
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VII. INFINITE -GAIN STATE-SPACE REALIZATIONS 



Although the idea of realizing high-order transfer functions by 
factoring them into first- and second-order product terms and realizing 
these using first- and second-order building-block networks has been 
mentioned in the previous chapters, it is in this chapter that the 
practicality and the advantages of such an idea are shown in detail. 

The realization methods in this chapter make use of state-space 
concepts to develop the basic second-order building block made up of 
integrators and summers, both of which use infinite-gain operational 
amplifiers. Consequently, such realization methods are classified 
under the category of Infinite-Gain Infinite-Gain Realizations. 



A. THEORETICAL DEVELOPMENT 

Consider the general open-circuit voltage transfer function 



V (s) 
o 

V^(s) 



a + a , s + 
o 1 



b + b, s + 
o 1 



^ n-1 ^ n 

+ a , s + a s 
n-1 n 

+ b , s + s 
n-1 



(7-1) 



A state representation of this transfer function is possible by breaking 
up equation (7-1) [Ref. 34] as follows: 



V (s) W (s) V (s) 

o _ F . o 

V^(s) Vj(s) Wp(s) 



where 
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1 



(7-3a) 



Wp(s) 

V^(s) 



b + b. s + 
o 1 



^ , n-1 ^ n 

. + b , s + s 
n-1 



and 



V (s) 
o 

Wp(s) 



n-1 n 

a +a,s + ,...+a ,s +as 
o 1 n-1 n 



(7-3b) 



where Wp(s) is the Laplace transform, with initial conditions set equal 
to zero, of a dummy variable, w(t) . Expanding equations (7-3a) and (7-3b) 
and getting the inverse Laplace transforms of the expanded equations 
give 



d w 



dt 



n 



= - b 



d*^ ^w 
n-1 ^^n-1 



, dw , . 

. - b^ — — b w + v^ 
1 dt o 1 



(7-4a) 



and 



^ dw ^ 
V = a w + a. ~ + 
o o 1 dt 



d w . d w 
. + a T — + a 

dt"-' " dt" 



(7-4b) 



where all variables are functions of time, 



If X, = w, x„ = dw , .... , and u = v, and y = v , two matrix 
' 2 — 1 o 



equations result: 



1 

1— 1 

1 




I 

O 


1 


0 


1 

O 




I—* 

1 




1 

o 

1 


€ 

^2 




0 


0 


1 ... 


0 




^2 




0 




= 




0 


• e • • 


• 




€ 


+ 


• 


X , 
n-1 




0 


0 


0 


1 




n-1 




0 


X 

n 




-b 

o 






1 

n-1 




X 

n 




1 


m 














_ 







(7-5a) 



H 
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and 




Equation (7-5a) may be written as 




X 

n- 



1 



X 

n 



(7-5b) 



X = A X + B u , (7-6) 

where x is the state vector, x is the derivative of the state vector, A 
is the matrix determined by the properties of the system and B u represents 
the inputs to the network. 

Equation (7-5b) may be expressed as 



y = C X + D u . (7-7) 

A signal flow graph can be drawn to represent equations (7-5a) and 
(7-5b) . This is shown in Figure 7-1. Each element of the graph can be 
realized by an integrator or a summer (inverter) , both of which are made 
up of operational amplifiers, resistors and capacitors, as shown in 
Figure 7-2 . 

If a high-order transfer function is to be realized by directly 

converting its signal flow graph as indicated above, the roots of such 
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Figure 7-1. Signal flow graph of general nth-order open-circuit voltage transfer function 
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Figure 7 -2a. A summer 




1S4 



high-degree polynomials can become very sensitive to the coefficients, 
so that very high accuracy is required in the resistive components that 
determine these coefficients. 

For this reason the most practical procedure is to factor the 
transfer function into first- and second-order product terms. For such 
terms, the high accuracy is not required and cascading the realizations 
of these product terms gives the transfer function in insensitive form. 



B. KERWIN, HUELSMAN, AND NEWCOMB [1967] 

Kerwin, Huelsman, and Newcomb [Ref. 35] recognized the value of 
the state-space approach and subsequently used its concepts to generate 
one of the first second-order building blocks . 

The general second-order transfer function in equation (2-14) is 
now considered: 



V. (s) 
V^(s) 



a + a, s + a.s 
o i z 

b + b, s + s^ 
o 1 



(7-8) 



If equation (7-8) is expanded as in Section A of this chapter using 
equation (7-2), the two resulting equations are: 



w = - b w - b,w + V- 
o 1 1 



V- = a w + a.w + a_w 
2 o 1 2 



(7-9a) 

(7-9b) 



If w = X, , w = X , v^ = u and v = y, Laplace transforming equations 

J. ^ J. 



(7-9) produces the following state equations: 
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(7-lOa) 



sXj(s) = X2(s) 

sX (s) = -b X (s) - b X (s) + U(s) (7-lOb) 

2 o 1 12 

Y(s) = a X (s) + a X (s) + a_sX (s) . (7-lOc) 

o i 12 2 2 

The signal flow graph corresponding to equations (7-10) is shown 
in Figure 7-3 . 

If each element of this second-order signal flow graph were re- 
placed by integrators and summers where appropriate, the resulting 
circuit would be the second -order building -block network of Kerwin, 
Huelsman and Newcomb. This is shown in Figure 7-4. 

Although four operational amplifiers, ten resistors and two capa- 
citors are required, this network can realize not only the general 

second -order transfer function in equation (7-8) by taking V as the 

0 

output voltage, but high-pass, bandpass and low-pass transfer functions 
as well, using as output V^, and respectively. 

By using the output equation of an integrator, which is 

V = - 1 V, . (7-lla) 

and that of a summer. 



V = 
o 




1 + R,G 
f n 



V . 
pi 



I 

^=1 



R.G V . 
f n ni 



(7-llb) 



with Figures 7-2a and 7-2b used as reference figures, where G and 

P 

G represent the sum of the input conductances of the positive and 
n 
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Figure 7-3. Signal flow graph of general second-order open-circuit voltage transfer function 
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Figure 7-4. Second-order building -block network of Kerwin, Huelsman, and Newcomb 



negative voltage inputs respectively, the different transfer functions 
for Figure 7-4 may be written in terms of the component elements as 
follows: 

General Second-Order Transfer Function 



RJR. + Rj 
0 o 

D(s) 



(7-12a) 



High-Pass Transfer Function 




D(s) 



Bandpass Transfer Function 



(7-12b) 




D(s) 



Low-pass Transfer Function 




K 



D(s) 



where 



^ 

Rj (R, + R^) 



/ 



(7-12c) 



(7-12d) 



(7-13) 
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and 



2 R, (R + R ) R C s + R 
D(s) = RgC^RgC^s . h 



(7-14) 



At the resonant frequency, Co , 

o 



CO = 
o 



R 






(7-15) 



for which R C R_C s + R/R = 0, 
o i y z o 



''3 r , 



R, 



(7-16) 



Use of Figure 7-4 and equation (7-12) is simplified if 






(7-17) 



such that the transfer functions in equation (7-12) become 




K, 



R7 (1 + V 



s + 1 + R^ s + 

1 + R„ 



s + 1+Rs+R 






(7-18a) 




K 



1 



2 

s 



2 

s 



+ 1 R 
1 + R, 



s + R 



(7-18b) 
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(7-18c) 





s 



2 

s 



+ 1 + R 
1 + R, 



s 



+ R 



V, 



K, 



s + 1 + R s + R 
1 + R„ 



where 






(IH-R^) 



In this normalized case, uj^ = \J~R and 



(7-18d) 



(7-19) 




0 




(7-20) 



At this point, an example illustrating the realization of a desired 
transfer function using the second-order building-block network in 
Figure 7-4 is appropriate. 

Example 7-1. Realize the bandpass filter transfer function of 
Example 2-4: 



H(s) = 



- 62.83 s 

s^ + 1.257 s + 3947.84 



(7-21) 



Since equation (7-21) was derived directly from the second-order 
bandpass transfer function form 



H(s) = 



- K CO s 
o 



2 2 
s + ( 4 (Jq/Q) s + 



(7-22) 
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the two variable element values can be determined from equations (7-17) 
and (7-20) as 

R = 3947.84 ohms (7-23a) 

^2 " "Tt fl" " • (7-23b) 



In the case where the realization of an odd-order transfer function 
is involved, the necessity for a first-order network becomes apparent. 
Such first-order sections can be realized to within a gain constant by 
lag-lead networks . 

The lag network shown in Figure 7-5a is described by the transfer 
function 



V R^Cs + 1 

O _ _2 

Vf “ (Rj+R2)Cs + 1 



(7-24) 



while the transfer function of the lead network shown in Figure 7-5b is 



V 



o 




'^2 

R1R2CS + Rj + R2 



C. HOLLENBECK [1969] 



(7-25) 



It is interesting to note that a simplified version of what is 
essentially the building-block network of Kerwin, Huelsman and Newcomb 
has already been put into production. Such a circuit, shown in Figure 
7-6, was discussed by Hollenbeck [Ref. 36]. 

The different transfer functions can be derived in a manner 
analogous to those in equations (7-12). They are as follows: 
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Figure 7-5a. First-order lag network 




Figure 7-5b. First-order lead network 
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Figure 7-6. A universal hybrid active filter 



For the high-pass output: 



V^(s) 



Fjd+F^) s‘ 



s + iUjFj (1 +F 2 ) s + co^ 



(7-2 6a) 



For the bandpass output: 






^ co^ CO/2 



(7-26b) 



For the low-pass output: 



^ ^ F, (1 ^F^) 



(7-26c) 



The parameters in equations (7-2 6) are defined in relation to circuit 
components as follows: 



^6 

F = ^ 

1 2R^ + R 



60 j “ 



1 



R7C7 



^2' 



R; 

r! 



= 



'^ 8'=8 






If the denominator coefficients of equations (7-22) and (7-26b) are 



compared, the following relations can be made: 






(7-28) 
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92.5 K{1 



S 10.81 
10 ^ 

R, = -^ = 1.125 kO. . (7-33) 

6 o 

D. SALERNO [1969] 

Using a slightly different approach, Salerno [Ref. 3 7] expanded 
the general voltage transfer function 



V 



o 




as follows: 



n n-1 

a s + a 1 s 
n n-1 



+ . 



u ri“l 
> + b , s + 
n-1 



+ a 
o 

+ b 

o 



(7-34) 



., / n , n-1 

V (s + b , s 
o n-1 



+ . . . + b ) = V, (a s^ 
o In 



+ a 



n-1 



n-1 

s 



+ . . . + a ) . (7-35) 

o 

By rearranging the terms of equation (7-35), the following general 
relation is arrived at: 



V = V, /a + a , 
o In n-1 



. + a 



o 

n 



- V /b - 
of n-1 



+ b 



jO 

n 



(7-36) 



The general second-order transfer function 



,, 2 

V a„s + a,s + a 

0 _ _2 1 o 

V ~ ~ 2 

1 s+b,s+b 

1 o 

can be similarly expressed as 



(7-37) 
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(7-38) 




Again using inverters whose gain is - 1 and integrators whose gain is 
- l/s, a schematic representation of equation (7-3 8) can be made, as 
shown in Figure 7-7. The circuit corresponding to this schematic is 
shown in Figure 7-8. The coefficients of the desired transfer function 
can be related directly to the components of the circuit by the following 
equations: 



a 

o 






a 



1 




a 



2 










(7-39) 



The other components are governed by the following relations: 




= l/C^ . (7-40) 

A similar procedure would realize a first-order transfer function 
if cascaded realizations of first- and second-order product terms are 
desired. An example of Salerno's procedure will clarify this point. 

Example 7-3. Realize the fifth-order Butterworth low-pass transfer 
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Figur6 7 7. Schorndtic rGpres8nt6tion of SdlGrno's r66liz3tion thG gonordl SGcond-ordGr tronsfor function 
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Figure 7-8. Salerno's circuit realization of the general second -order transfer function 



function of Example 2-2: 



V 

_c 

V, 



1+s 

. (Jo 



1+0. 618s + ^ ^ ^ 



, (7-41) 



CJo to n 

o 



1+ 1 ,618s + ^ 

U) u 

o o 



where = 2 TV (10). 
^o 



Equation (7-41) can be expressed as three cascaded transfer 



functions: 



V. 



V, 



- 1 



1 + s 



U) 



(7-42a) 



V, 

V. 



1 + 0 . 618s + s 



Ld, 



(d. 



(7-42b) 



V 

_c 

V, 



1 + 1. 618s + s 



U) 



to. 



(7-42c) 



Expanding equation (7-42a) gives 



''2 = 



00 o - to o 

s s 



(7-43) 



the schematic and circuit realization of which are shown in Figure 7-9. 



V = OOo \j - \j f CJq + 0«618 



2 2 3^3 



(7-44) 



cO 



V = ^ o V„ - V 

o ;; — 3 o 



1.618 Co 

^ o + o 



(7-45) 



171 



COo 




— J 

Figure 7-9a. Schematic representation of a first-order transfer 
function 
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Figure 7-10. Schematic representation of a second-order low-pass 
transfer function 







Figure 7-11. Circuit realization of a second-order low-pass 
transfer function 
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Figure 7-12 . Circuit realization of Example 



Xj = - ax 



j - kj / b ' I ma - c 



V, 



(7-48a) 



S ^ b X 
k ^ 



(7-48b) 



V 



J_ sgn (ma - c) x + d _ ^x ^4-n,v^,(7-48c) 



1 k^ |ma - c| 



where k and k are arbitrary positive real numbers and where 

X ^ 



sgn (K) = + if K is positive 

= - if K is negative . (7-49) 

In relation to the notation in Ref. 34, x^^ = w and x^ = w . The circuit 
configuration realizing equation (7-47) and the set of state equations in 
equation (7-48) is shown in Figure 7-13. The element -coefficient 
relations for this configuration are as follows: 



R, = 1 



aC, 






R, 



fbC. 



^6 ^3 



= 









ma - c C, 



*^7 P 2 '^10 



«8 = 



|ma - c \fT R 
[mb - dj 



10 



■^9 = 



m 



R^q for Cases A, B, C 



= b 



Rjq for Case D 



(7-50 
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Figure 7-13. Tow's second-order building-block network 
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Figure 7-14. Tow's active bandpass network 
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Figure 7-1 5. Tow's active low-pass network 



To minimize change in pole locations due to finite gain of the 
operational amplifiers, k and k are generally chosen to be unity but 

X 

in cases where it is desired to minimize the spread in element values, 
other convenient values may be chosen. 

For cascaded realizations, Tow also specified a method for realizing 
the first-order product term of the form 

^ ^ K . (7-54) 

Vj s + a 

Direct coefficient-matching of equation (7-54) with 



V 



o 




1 

RC 

s + R + R, 
in 

RR. 




(7-55) 



in 

gives the values of R and C, where R, C, and R.^ are as indicated in 
Figure 7-16 . 

Example 7-4. Realize the fourth-order bandpass transfer function 
of Example 2-4: 



H(s) = 



1.18 s 

(s^ + 0.94 s + 3913.80) (s^ + 0.95s + 3982.18) 



(7-56) 



For a cascaded realization, equation (7-56) can be factored into 



two product terms as follows; 



H(s) - 



1.18s A 

s^+0.94s + 3913. 80y 



. (7-57) 



s + 0.95s + 3982.18 
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Figure 7-16. Cascaded first- and second-order sections where dotted block indicates 
realization of first-order pole 



The first product term is first realized using the coefficient-element 



relatioiiS in equation (7-50) and noting that the required second-order 
section is bandpass and hence, does not require the summing amplifier, 
A , as in Figure 7-14. Arbitrarily choosing C = C = 1 F and 

Hi • X ^ 

k = k = 1, the other element values are found to be: 

J. ^ 



''l 


= 10« = 


1.063 mH 


0.94 






= 10® 


15.99 k £1 


L 


|3913.80‘ 






= 10« = 


15.99 kSI 


3 


/3913.80' 






= 10® = 


847.0 Kll 




1.18 






= ^ = 


15.99 kO. 



(7-58) 



Cimilarly, for the same choice of , C„ , k , and k , the element 

X ^ X ^ 

values of the second-product term section are calculated to be: 



''i 


= 10® = 


1.052 


Mfl 


0.95 








= 10® = 


15.86 


kH 


2 


/3982.18 








= 10® 


15.86 


kQ 


3 


/3982.18 






'^4 


= 10® = 


1 mQ. 




''s 


= '^6 = "3 = 


15.86 


kH 
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F . ADVANTAGES AND DISADVANTAGES 

The state-space approach used in the filter transfer function 
realizations in this chapter has the distinct advantage of using well- 
developed theories from system and control applications. It takes 
advantage of direct analog computer simulation. In addition, the concept 
of using first- and second-order building-block networks significantly 
increases the stability and insensitivity of such realizations to fluctua- 
tions in the values of the passive RC components and in the gains of 
the operational amplifiers used. These sensitivity considerations are 
well discussed in Refs. 3 5 and 40. The building-block networks them- 
selves are versatile (low-pass, bandpass and high-pass outputs) and 
readily adaptable to cascade realization because of their low-output 
impedance. And most important from economic considerations, such 
building-block networks are compatible to the microminiaturization and 
mass production concepts of present-day technology. 
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VIII. CONCLUSION 



The synthesis techniques presented in the last four chapters are 
by no means the only ones that are possible with the use of operational 
amplifiers. It has been shown [Ref. 41] that the operational amplifier 
can be used to realize the three other classes of active devices - 
controlled sources, the negative-immittance converter and the gyrator - 
and consequently all the realization procedures which are applicable 
using these active devices are possible using operational amplifiers. 
Reference 42 gives an excellent treatment of such realization procedures . 
Operational amplifiers may even be used to take advantage of highly- 
developed passive filter synthesis procedures by direct inductance 
simulation [Refs. 43 - 45]. 

The use of operational amplifiers, therefore, offers a filter 
designer many advantages, among which is a wide choice of realization 
techniques. The particular method chosen is, of course, still highly 
dependent on the specifications of the particular filter desired. It has 
been pointed out, however, that the realization techniques using the 
state-space approach in Chapter VII of this paper offer much in the way 
of doing away with the expensive process of custom-building filters. 

One limitation of all realization techniques using operational 

amplifiers should be stressed at this time. Operational amplifiers 

typically have gain-bandwidth products of the order 1 MHz. Consequently, 

at frequencies much above a few kilohertz, the infinite-gain approximations, 
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on which feedback synthesis is based, break down. It should be 
emphasized, then, that realization techniques using operational 
amplifiers are limited to frequencies within and below the audio range. 
However, the fact that inductors become unwieldy at these very same 
frequencies brought the development of active RC filters in the first 
place . 

In conclusion, it is hoped that this paper has served to provide 
an orderly transition from the concepts of modern filter theory to a 
better understanding of the synthesis techniques involved in the real- 
ization of active RC filters using operational amplifiers. 
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APPENDIX A 



TOW’S FREQUENCY TRANSFORMATION PROCEDURE 

Reference 9 contains the normalized low-pass poles and zeros up 
to n = 9 of four different characteristic types - Butterworth, Chebyshev, 
Inverted Chebyshev and Cauer Parameter. Tow's transformation pro- 
cedure converts these poles and zeros into the poles and zeros of a 
desired bandpass or band-reject transfer function. 

Let the normalized low-pass complex pole or zero from Ref. 9 be 



no 



(Tn*' 



(A-1) 



For a desired bandpass or band-reject filter, either the resonant 

or center frequency, f , or the lower and upper frequency limits, f, and 

o 1 

f^, are normally specified. The following relations can consequently be 
determined: 



= 2TTf^ or 2TT(^) 



BW = f - £, 
u 1 



X = BW/f = 1/Q 
o 



(A-2) 

(A-3) 

(A-4) 



Equations (A-1) to (A-4) may now be used in the following low-pass- 
to-bandpass root transformation: 



•k k 

(s “ Sj)(s - s^ )(s - S 2 )(s - s^ ) (A-5) 
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where s is the complex frequency variable of the frequency-transformed 
transfer function and where 



Sj = - 1/2 ( cr + v) CaJ^ + j 1/2 {otJ + u) 



S 2 = - 1/2 ( cT - v) + j 1/2 ( o; - u) cO^ 



(A-6) 

(A-7) 



2 ^ 2 , v2 

K = 4 rr f - fj 

u 1 



(A-8) 



and 



if = 



(A-9) 



X 

n 

= \l/2 {4 - 



" -t 



(A- 10) 



+ ^4 - 






1/2 



(A- 11) 



V - CT U )/ u 



(A- 12) 



* "k 

s and s are the complex conjugates of s and s respectively. The 

X ^ X ^ 

normalized -low-pas s-to-band -reject root transformation is: 



(S + ^ - } u> ) is + ^ + j cj) 
n*^n nn^^n n 



•k * 

(s - Sj)(s - s^ )(s - S^){s - ) 

/ 2 ^ TTi 
(s + CJ ) 

o 



(A- 13) 



where s , s , u, and v are as given by equations (A-6), (A-7), (A-11) 

X ^ 



and (A-12) respectively and 
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= 



(A- 14) 



U) = 



^ n 

, - 2 . 2>2 

' ^ > 



OJ ^ 

n 

^ 2 2,2 



(A- 15) 



189 



APPENDIX B 



COMPUTER PROGRAM 1 



SUBROUTINE FREQR ( NCN,NCD, A, B , FMI N , FM AX , X , Y I 



PURPOSE - TO PLOT THE FREQUENCY RESPONSE OF ANY GIVEN 
TRANSFER FUNCTION WITHOUT POLES ON THE JW-AXIS, 

NCN - THE NUMBER OF NUMERATOR COEF FI C I ENTS-, 

NCD - THE NUMBER OF DENOMINATOR COEE E I Cl ENTS. 

A - THE ARRAY OF NUMERATOR COEFFICIENTS STARTING FROM THE 
THE LOWEST-DEGREE tfrMo 

B - THE ARRAY OF DENOMINATOR COEFFICIENTS STARTING FROM THE 
THE LOWEST-DEGREE TERM. 

FMIN - THE LOWER FREQUENCY LIMIT IN HERTZ TO BE PLOTTED- 
FMAX - THE UPPER FREQUENCY LIMIT IN HERTZ TO BE PLOTTED, 

X - THE ARRAY OF FREQUENCY POINTS IN HERTZ TQ PE PLOTTED. 

Y - THE ARRAY OF MAGNITUDE POINTS TO BE PLOTTED. 

ADDITIONAL SUBROUTINE REQUIRED 



SUBROUTINE DRAW FROM THE NPS WoRo CHURCH COMPUTER 
CENTER LIBRARY 



SUBROUTINE FREQR ( NCN , NCD , A , B , FMI N , FMAX , X , Y ) 

DIMENSION A(30 ) , B( 30) , X( 201) ,Y( 201 I 
COMPLEX H, StD, BI 

REAL=!'8 ITITLE(12)/» FIRST LINE OF DESIRED TITLE, 
o SECOND LINE OF DESIRED TITLE fqp PLOT.*/ 

REAL LABEL/4H / 

WRITE! 6, 105) 

WRITE(6,1C4) ( A( J) , J=1 ,NCN) 

WRITE(6,103) 

WRITE(6,104) (B( J) ,J = 1 ,NCD) 

WMIN=FMIN46. 283184 
WMAX=FMAX*6o 2831 84 
DELW=( WMAX-WMIN) /200. 

DO 1000 M=l,201 
W=WMIN+DELW4(M-1 ) 

S=CMPLX(OoO,W) 

H=(0o0,0.0 ) 

D=(0o0,0o0) 

DO 1 1=1, NCN 
AI=A(NCN+1-I ) 

BI=CMPLX(AI ,0o0) 

H=H*S+BI 

1 CONTINUE 

DO 2 1=1, NCD 
AI=B(NCD+1-I ) 

BI=CMPLX(AI ,OoO) 

D=D«S+BI 

2 CONTINUE 
H=H/D 

X(M)=AIMAG(S)/6o 283184 
X( M) =X( M)-WM IN/6 o 2 83184 
Y(M)=CABS(H) 

1T)C0 CONTINUE 

CALL DRAW! 2CC , X , Y, 0 , 0 , L AB E L , ITITLE,C,0,Q,0,C,C,5,8,i, 
♦ LAST ) 

105 FORMAT! IHl , lOX, • THE COEFFICIENTS OF THE NUMERATOR ARE’ 

*,///) 

103 .FORMAT! ///, lOX, * THE COEFFICIENTS OF THE DENOMINATOR 
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104 



♦ ARE* »/// ) ^ 

F0RMAT(10XtF20o5,5X,F20o5 
RETURN 
END 



/ ) 
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APPENDIX C 



NATHAN'S METHOD OF MATRIX ANALYSIS 

Consider the n-terminal network with an operational amplifier 
connected between its i th and j th terminals as shown below: 




If the network is considered by itself, without the operational 
amplifier, it can be described by the following matrix equation: 

[I] = [Y] [E] , (C-1) 

where [I] is the column matrix of terminal currents, [E] is the column 
matrix of terminal voltages, and [Y] is the admittance matrix for the n- 
terminal network. 

If the network is now considered taking the presence of the opera- 
tional amplifier into consideration, several simplifying operations can be 
performed . 

Since E. = 0 by the concept of virtual earth, E^ may be removed 
from [E] , cancelling effectively the i th column of [Y] . 

Also, since L will always be of such value to satisfy the con- 
straint that E. = 0, I. may be eliminated, again effectively cancelling 
1 3 

the j th row of [Y] . 
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These operations yield a new matrix equation: 



[r] = [Y’] [E’] , (C-2) 

where [!'] is [I] with I. deleted, [E‘] is [E] with E. deleted, and [Y'] is 

] 1 

[Y] with the j th row and the i th column deleted . 

If the voltage transfer function between two terminals, p and q, of 
the n-termi*ial network is desired, matrix algebra yields: 

[E’] = [Y-] [r] . (C-3) 

The voltage at terminal p can be written as 

n 

E = 1 > Y^*^ I , (C-4) 

" |Y’I ^ 

where the Y^"^ terms are cofactors associated with corresponding terms 
in [Y’]. 

If all terminals except p, i, and j are left open-circuited and if 

terminal p is excited by a voltage source, the only non-zero terminal 

currents will be I and I.. I. = 0 because of the operational amplifier's 
P 3 1 

infinite input impedance . Application of these conditions into equation 
(C-4) gi'. es 

E = Y^^ I . (C-5) 

P |y'| p 

Under these same conditions. 



E 

q 




(C-6) 
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Therefore, provided p is not driven and q is not driving. 



E 

_a 

E 

P 



Y 



pq 



yPP 



.(C-7) 
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